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Abstract. This paper develops and implements a new algorithm for calculating wave 
trace invariants of a bounded plane domain around a periodic billiard orbit. The algorithm 
is based on a new expression for the localized wave trace as a special multiple oscillatory 
integral over the boundary, and on a Feynman diagrammatic analysis of the stationary phase 
expansion of the oscillatory integral. The algorithm is particularly effective for Euclidean 
plane domains possessing a Z2 symmetry which reverses the orientation of a bouncing ball 
orbit. It is also very effective for domains with dihedral symmetries. For simply connected 
analytic Euclidean plane domains in either symmetry class, we prove that the domain is 
determined within the class by either its Dirichlet or Neumann spectrum. This improves 
and generalizes the best prior inverse result (cf. [Zl[|Z21llSZp that simply connected analytic 
plane domains with two symmetries are spectrally determined within that class. 



1. Introduction 

This paper is part of a series (cf. |Z5l IZ4] ) devoted to the inverse spectral problem for 
simply connected analytic Euclidean plane domains Q. The motivating problem is whether 
generic analytic Euclidean drumheads are determined by their spectra. All known counterex- 
amples to the question, 'can you hear the shape of a drum?', are plane domains with corners 
[GWWlj . so it is possible, according to current knowledge, that analytic drumheads are spec- 
trally determined. Our main results give the strongest evidence to date for this conjecture 
by proving it for two classes of analytic drumheads: (i) those with an up/down symmetry, 
and (ii) those with a dihedral symmetry. This improves and generalize the best prior results 
that simply connected analytic domains with the symmetries of an ellipse and a bouncing 
ball orbit of prescribed length L are spectrally determined within this class \Z1[ IZ21 llSZj . 

The proofs of the inverse results involve three new ingredients. The first is a simple and 
precise expression (cf. Theorem 13.11) for the localized trace of the wave group (or dually the 
resolvent), up to a given order of singularity, as a finite sum of special oscillatory integrals 
over the boundary dQ of the domain with transparent dependence on the boundary defining 
function. Theorem 13.11 is a general result combining the Balian-Bloch approach to the wave 
trace expansion of [Z5] with a reduction to boundary integral operators explained in [Z4] . 
Presumably it could be obtained by other methods, such as the monodromy operator method 
of lantchenko, Sjostrand and Zworski [SZt IISZ|| . Aside from this initial step, this paper is 
self-contained. 

The next and most substantial ingredient is a stationary phase analysis of the special 
oscillatory integrals in Theorem 13.11 To bring order into the profusion of terms in the wave 
trace (or resolvent trace) expansion, we use a Feynman diagrammatic method to enumerate 
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the terms in the expansion. Diagrammatic analyses have been previously used in [AG] (see 
also [BuJ) to compute the sub-principal wave invariant. A novel aspect of the diagrammatic 
analysis in this paper is its focus on the diagrams whose amplitudes involve the maximum 
number of derivatives of the boundary in a given order of wave invariant. A key result, 
Theorem 14.21 is that only one term, the principal term in Theorem 13. contributes such 
highest derivative terms. That is, the stationary phase expansion of the principal term 
generates all terms of the jth order wave invariant (for all j) which depend on the maximal 
number 2j — 2 of derivatives of the curvature of the boundary at the reflection points. In the 
principal term, the 'transparent dependence' of the phase and amplitude on the boundary 
is encapsulated in the simple properties of the phase and amplitude stated in the display 
in Theorem 14.21 Only these properties are used to make the key calculations of the wave 
invariants stated in Theorem 15.11 

This focus on highest derivative terms in each wave invariant turns out to be crucial for 
the inverse spectral problem on domains with the symmetries studied in this article. The 
third key ingredient is the analysis in ^ of these highest order derivative terms in the case 
of domains in our two symmetry classes. The main result is that the other terms in the wave 
invariants are redundant, and further that the domain can be determined from the wave 
invariants within these symmetry classes. These results are based on the use of the finite 
Fourier transform to diagonalize the Hessian matrix of the length function, and an analysis 
of Hessian power sums. 

As this outline suggests, we take a direct approach to calculating wave trace invariants 
and do not employ Birkhoff normal forms as in [Gl [ZTl [Z2l [Z3l IISZ] . We do this because 
the classical normal form of the first return map does not contain sufficient information 
to determine domains with only one symmetry. Therefore one would need to use the full 
quantum Birkhoff normal form. But we found the calculations based on the Balian-Bloch 
approach simpler than those involved in the full quantum Birkhoff normal form. 

1.1. Statement of results. Let us now state the results more precisely. We recall that the 
inverse spectral problem for plane domains is to determine a domain Q as much as possible 
from the spectrum of its Euclidean Laplacian in fl with boundary conditions B: 



The boundary conditions could be either Dirichlet B(p = (p\dn, or Neumann B(p = di^iplg^ 
where di, is the interior unit normal. 

We briefly introduce some other notation and terminology, referring to §2 and to [KTJ 
-(PS] for further background and definitions regarding billiards. By Lsp{Q) we denote the 
length spectrum of Q, i.e. the set of lengths of closed trajectories of its billiard fiow. By 
a bouncing ball orbit 7 is meant a 2-link periodic trajectory of the billiard fiow. The orbit 
7 is a curve in S*Q which projects to an 'extremal diameter' under the natural projection 
vr : S*Q — s> Q, i.e. a line segment in the interior of Q which intersects dQ orthogonally at 
both boundary points. For simplicity of notation, we often refer to vr(7) itself as a bouncing 
ball orbit and denote it as well by 7. By rotating and translating Q we may assume that 
7 is vertical, with endpoints at A = (0,^) and B = (0, — -1). In a strip T^{AB) of width 



(1) 
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y = f-{x) 
(0,-f) 

Figure 1 . as a pair of local graphs 

epsilon around 7, we may locally express dVt = dQ^ U dQ~ as the union of two graphs over 
the X-axis, namely 

(2) dn+ = {y = f+{x), xe{-e,e)}, dQ' ^ {y = f.{x), x e {-e,e)}. 

Our inverse results pertain to the following two classes of drumheads: (i) the class Vi^l of 
drumheads with one symmetry a and a bouncing ball orbit of length 2L which is reversed 

by <t; and (ii) the class Pm,L (^^ ^ 2) of drumheads with the dihedral symmetry group 
and an invariant m-link reflecting ray. Let us define the classes more precisely and state the 
results. 

1.1.1. Domains with one symmetry. The class consists of simply connected real-analytic 
plane domains Q satisfying: 

• (i) There exists an isometric involution cr of Q which 'reverses' a non-degenerate 
bouncing ball orbit 7 — > of length — 2L. Hence f+{x) — — 

• (ii) The lengths 2rL of all iterates 7'' (r = 1, 2, 3, . . . ) have multiplicity one in Lsp{fl), 
and in the elliptic case, the eigenvalues e*" of the linear Poincare map satisfy that 
a — —2 cos I does not belong to the 'bad set' i3 = {a = 0, —1, 2, —2}. 

• (iv) The endpoints of 7 are not vertices of dQ. 

Let SpecB(Jl) denote the spectrum of the Laplacian Aq of the domain Q with boundary 
conditions B (Dirichlet or Neumann). 

Theorem 1.1. For Dirichlet (or Neumann) boundary conditions B, the map SpecB : T^i,l ^ 
is 1-1. 

Let us clarify the assumptions and consider related problems on Z2-symmetric domains: 

(a) Under the up-down symmetry assumption, f+{x) — —f-{x) (see Figure (2)). Hence there 
is 'only one' analytic function / to determine. It is quite a different problem if a preserves 
orientation of 7 (i.e. flips the domain left-right rather than up-down), which amounts to 
saying that f± are even functions but does not give a simple relation between them. 

(b) Condition (ii) on the multiplicity of 2L means that 7 is the only closed biUiard orbit 
of length 2L. Since 7 = 7"^ for a bouncing ball orbit, the multiplicity is one rather than 
two. The method we use to calculate the trace combines the interior and exterior problems, 
and so one might think it necessary to assume that no exterior closed billiard trajectory (in 




STEVE ZELDITCH 



Figure 2. A domain in Vi^l 



the complement VL'^ of Vt) has length 2L. However, it is known that there exists a purely 
interior wave trace (cf. §1.2p and that the wave trace invariants at 7 are spectral invariants; 
we use the interior /exterior combination only to simplify the calculation. Therefore, it is 
not necessary to exclude exterior closed orbits of length L. When making stationary phase 
calculations, we only consider the interior closed orbits. 

(c) The linear Poincare map is defined in ^ In the elliptic case, its eigenvalues {e^*"} are 
of modulus one and we require that a = — 2 cos ^ lies outside the bad set B. In the hyperbolic 
case, its eigenvalues {e^"} are real and they are never roots of unity in the non-degenerate 
case. These are generic conditions in the class of analytic domains. We refer to the angles 
a as Floquet angles. The set B consists of angle parameters where certain functions fail to 
be independent as one 'iterates' the geodesic 7. The role of this set will be described more 
precisely in §1.2.31 

(d) Assumption (iii) is equivalent to f±\Q) 7^ 0. The third derivatives /]f'*(0) of f± at 
the endpoints of the bouncing ball orbit appear as coefficients of certain terms in the wave 
invariants, and we make assumption (iv) to ensure that the corresponding term does not 
vanish. Geometrically, f±\^) = only if the endpoints of the bouncing ball orbit are 
vertices of dVL, i.e. critical points of the curvature. This is a technical condition which we 
believe can be removed by an extension of the argument, as will be discussed at the end of 
the proof. We do not give a complete argument for the sake of brevity. 

As a corollary, we of course have the main result of [Zlt IZ21 IISZ] that a simply connected 
analytic domain with the symmetries of an ellipse and with one axis of a prescribed length 
L is spectrally determined within this class. 

Corollary 1.2. Let T>2 he the class of analytic convex domains with central symmetry, i.e. 
the symmetries of an ellipse. Assume that {rL^} are of multiplicity one in Lsp{Q) up to 
time reversal (r = 1, 2, 3, . . Then Specs: T> M.^ is 1-1. 

We give a new proof at the start of ^ since it is much simpler than the one-symmetry 
case and since the proof is simpler than the ones in \L1\ IZ2j . 
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y = f{x) 




Figure 3. Z2 symmetric domain with comers 

This inverse result is also true for non-convex simply connected analytic domains with 
the symmetries of the ellipse if we assume one axis has length L and is of multiplicity one. 
We stated the result only for convex domains because, by a recent result of M. Ghomi [Ghj . 
the shortest closed trajectory of a centrally-symmetric convex domain is automatically a 
bouncing ball orbit, hence it is not necessary to mark the length L of an invariant bouncing 
ball orbit. 

Theorem (II. ip removes the (left /right) symmetry from the conditions on the domains 
considered in \L1\ IZ2j . The situation for analytic plane domains is now quite analogous to 
that for analytic surfaces of revolution [Z3j, where the rotational symmetry implies that the 
profile curve is up/down symmetric but not necessarily left/right symmetric. 

Theorem 11.11 admits a generalization to the special piecewise analytic mirror symmetric 
domains with corners which are formed by reflecting the graph of an analytic function y = 
f{x) around the x-axis. More precisely, let f{x) be an analytic function on an interval [—a, a] 
(for some a) such that /(a) = f{—a) = and that / has no other zeros in [—a, a]. Then 
consider the domain Qf bounded by the union of the graphs y = ±f{x). 

Let JF be the class of real analytic functions with the stated properties, and consider those 
/ for which precisely one critical value of / equals L/2. The vertical line through [x, ±L/2) 
is then a bouncing ball orbit. We further impose the same generic conditions on Qf as in 
Theorem 11.11 We denote the resulting class of real analytic graphs by J-'l- 

Theorem 1.3. Up to translation (i.e. choice of a), the Dirichlet (or Neumann) spectrum 
ofQf determines f within J-'l, i.e.: Spec: Tl ^ i-i- 

The proof is identical to that of Theorem 1 1 . 1 1 once it is established that there exists a wave 
trace expansion around the length t = 2L of the bouncing ball orbit for domains in T with 
the same coefficients as in the smooth case. This fact follows from work of A. Vasy [Vj on 
the Poisson relation for manifolds with corners. In other words, the presence of corners does 
not affect the wave trace expansion at the bouncing ball orbit. 

1.1.2. Dihedrally symmetric domains. The second class of domains is the class T>m,L of di- 
hedrally symmetric analytic drumheads Q, i.e. domains satisfying: 

• (i) tQ = Q for all r G Dm] 

• (ii) Dm leaves invariant at least one m-link periodic reflecting ray 7 of length 2L; 

• (iii) The lengths 2rL have multiplicity one in Lsp{Q) 
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Figure 4. A Z^s-symmetric domain 

We then have: 

Theorem 1.4. For any m>2, Specs ■ T^m,L ^ IR+ is 1-1- 

We recall that Dm is the group generated by elements {a, -R27r/m} where -R27r/m is counter- 
clockwise rotation through the angle 27r/m and where cx^ = 1, with the relations aR2n/nO' = 
R -2n/n- Also, by an m-link periodic reflecting ray we mean a periodic billiard trajectory with 
m points of transversal reflection off dQ. It is easy to see that such a ray exists if Q is convex. 
In general, it is a non-trivial additional assumption. With this proviso. Theorem (11. 4p is a 
second kind of generalization of the inverse spectral result of \Z1\, |Z2 J for the class P2,l of 
'bi-axisymmetric domains'. That result obviously covers the classes T>2n,L, but the general 
case is new. For any prime p, the result for T'p,L is independent of any other case where p 
does not divide n. 

1.2. Overview. Let us give a brief overview of the proofs. 
We denote by 

E2it,x,y) = ^costXj(pj{x)(pj{y) 
j 

the kernel of the even part of the wave group cost^/A^, generated by the Laplacian A§ 
of ([T]) with either Dirichlet Bu = u\gn or Neumann Bu = d^ulgn boundary conditions. Its 
distribution trace is defined by 

„ oo 

(3) TrlnE^it) := / x, x)c/x = V costAj 

Jn 

When is the length of a non-degenerate periodic reflecting ray 7 of the generalized 
billiard flow, and when the only periodic orbits of length are 7 and 7"^ (the time- 
reversal of 7), then TrlQE^it) is a Lagrangian distribution in the interval {L^ — e, L^ + e) for 
sufficiently small e, and has the following expansion in terms of homogeneous singularities: 
(see |GM] . Theorem 1, and also page 228; see also [PS\ Theorem 6.3.1). 

Let 'J be a non- degenerate billiard trajectory whose length is isolated and of multiplicity 
one in Lsp{Q). Then for t near L^, the trace of the even part of the wave group has the 
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singularity expansion 
(4) 

oo 

TrluE^it) ~ ^{a^{t-L^+iO)-^ + a^o^og{t-L^ + iO) + ^a^k{t-L^ + iO)Hog{t-L^ + iO)}, 

k=l 

where the coefficients a^^ (the wave trace invariants) are calculated by the stationary phase 
method from a microlocal parametrix for at 7. 

Here, is a sum of the contributions from 7 and 7"^, which are the same. In general, 
the contribution at t = is the sum over all periodic orbits of length L^. The sum to the 

right of 3fJ is the trace of the wave group e**v^^; the trace of the even part E^it) of the wave 
group equals the real part of that trace. 

In |Z5] . §3.1, this expansion was reformulated in terms of a regularized trace of the interior 
resolvent R^{k + iT) = -(Ag + (fc + ir)^)-! : H'{n) H'+^{n), with A; G M, r > and with 
boundary condition B. The Schwartz kernel or Green's kernel G^{k + iT,x,y) G V{Q x Q) 
of the resolvent is the unique solution of the boundary problem: 

r -(Ag + {k + iTf)G'^{k + IT, X, y) = 6y{x), (x, yen) 

(5) 

[ BG%{k + iT, x,y)=0, x G dQ 

Let p G G^{L^ — e, + e) be a cutoff, equal to one on an interval (L^ — e/2, + e/2) 
which contains no other lengths in Lsp(fi) occur in its support, and define the smoothed 
(and localized) resolvent with a choice of boundary conditions by 

(6) R^^{k + iT):= / p{k - p){fi + iT)R^{fi + iT)dfi. 
The definition is chosen so that 

roo 

(7) RUk + tT)= / p(t)e^('=+-)*E^(t)rft. 

Jo 

Then the smoothed resolvent trace admits an asymptotic expansion of the form 

00 

(8) TrlnRepik + ir) ~ VB,^{k + ir) ^ B^jk-^, k ^ 00, 

j=0 

where 

• Vs^-yik + ir) is the symplectic pre-f actor 



VB,-y{k + iT) = Go £5(7)- 



V|det(/-P,)| 

is the Poincare map associated to 7 (see ^for background); 
£5(7) is the signed number of intersections of 7 with dQ (the sign depends on the 
boundary conditions; ±1 for each bounce for Neumann/Dirichlet boundary condi- 
tions) ; 

is the Maslov index of 7; 
Gq is a universal constant (e.g. factors of 2tt) which it is not necessary to know for 
the proof of Theorem 11.11 
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The resolvent trace (or Balian-Bloch) coefficient B^j associated to a periodic orbits 7,7"^ is 
easily related to the wave trace coefficient a^^k- We henceforth work solely with the expansion 
([8]), which we term the 'Balian-Bloch expansion' after [BB2|. In fact, we actually analyze 
the closely related resolvent trace asymptotics along logarithmic curves k + ir log k in the 
upper half plane. It is clear that the 'Balian-Bloch coefficients' B^ j are spectral invariants 
and it is these invariants we use in our inverse spectral results. 

As mentioned above, the inverse results have three main ingredients, which we now describe 
in detail as a guide to the paper and its connections to [Z4l IZ5j . 



1.2.1. Reduction to boundary oscillatory integrals of the wave trace. The ffist step (Theorem 
13. ip is a reduction to the boundary of the wave trace. This reduction was largely achieved in 
\Z5\ IZ4] by means of a rigorous version of the Balian-Bloch approach to the Poisson relation 
between spectrum and closed billiard orbits |BBlt[BB2] . It expresses the wave trace localized 
at the length of a periodic reflecting ray, up to a given order of singularity, as a finite sum 
of oscillatory integrals I'^^p{k + ir) over the boundary (see f|T9l) . It is related in spirit to the 
monodromy operator approach of |SZl IISZ] . 



1.2.2. Feynman diagram analysis and proof of Theorem \4-^ The second ingredient is a 
stationary phase analysis of the oscillatory integral expressions for the wave invariants at 
transversally reflecting periodic orbits. The key role is played by a (Feynman) diagrammatic 
analysis of the stationary phase expansions, which has not previously been used in inverse 
spectral theory (see [AG] for prior use in calculated the sub-principal invariant). As reviewed 
in §5.11 the terms of stationary phase expansion correspond to labelled graphs F and the 
coefficients of the stationary phase expansion can be expressed as 'Feynman amplitudes' 
determined by the graphs F. The Euler characteristic of F corresponds to the power k~^ of 
k in the wave trace expansion. 

The inverse spectral problem involves a novel point of the diagrammatic analysis: namely, 
to separate out the (labelled graphs) of Euler characteristic —j whose amplitudes contain 
the maximum numbers (2j + 2, 2j — 1) of derivatives of dVL. In Theorem 14. 21 we prove that the 
terms in a given wave invariant which contain the maximal number of derivatives of dVL only 
arise in the stationary phase expansion of one principal term and its time reversal, whose 
amplitudes have special properties stated in table in Theorem 14.21 The principal terms are 
defined in Definition 14. 3[ Only the special properties of the phase and amplitude are used 
in the calculation of the wave trace invariants. 

The analysis leads to the explicit formulae for the top derivative parts of the wave in- 
variants at iterates of bouncing ball orbits in Theorem 15. 1[ For instance, in the symmetric 
bouncing ball case there is only one important diagram for the even derivatives /'•^-'■'(O) and 
two important diagrams for the odd derivatives /'•^-'^^''(O). Modulo terms involving < 2j — 2 
derivatives, the wave trace (or more precisely resolvent trace) invariants (cf. (jl])-® B^j-i 
take the form (cf. Corollary 15.111) : 

= {ALr)Arm^-'{2{wg2,,o) {hllrf^^^KO) 
(9) +4(^,.-.,3,o) (/.i^)^-,3^(/(^)(0)/(^-^)(0)) 

+4{w^. ..o) {hiiy-'ElUhltm^'Ko)f('^-'\o))}. 
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Here and throughout the paper we use the following notational conventions: 

• matrix elements of the inverse of the Hessian of the length function 
£ in Cartesian graph coordinates at Y' (cf. 

• Ar{0) is an f2-independent (non-zero) constant obtained from amplitude of the prin- 
cipal terms at the critical bouncing ball orbit. 

• Wq2],o (etc.) are certain non-zero combinatorial constants associated to Feynman 

graphs denoted here by etc. For a given graph Q, wg = where |ylut(^)| 

is the order of the symmetry group of the graph; see the discussion after (1561) . 

The amplitude value Ar{0) and the Wick constants may be evaluated explicitly. However 
it is not necessary for the proof of Theorem 11.11 to do so and it seems more illuminating to 
specify the origins, rather than their values, of the various constants. We note that the h^^. 
depend on, and only on, r and the eigenvalues of the Poincare map (i.e. on the Floquet 
angles) and on the length of 7. We also note that 7 = 7"^ when 7 is a bouncing-ball orbit 
(such an orbit is called reciprocal). 

The analysis shows that the non-principal oscillatory integrals only give rise to sub- 
maximal derivative terms in the wave invariants, completing the proof of Theorem 14.21 

1.2.3. Inverse results. The third ingredient is the analysis of the top derivative terms in the 
wave trace invariants in the symmetry classes above. The key point is determine the 2j — 1st 
and 2jth Taylor coefficients of the curvature at each reflection point from the j — 1st wave 
trace invariant for 7 and its iterates 7^. 

We note that the previously known inverse result for analytic domains with the symmetry 
of an ellipse drops out immediately from (Q, since the odd Taylor coefficients are zero. On 
the other hand, there is an obstruction to recovering the Taylor coefficients of / when there 
is only one symmetry: namely, we must recover two Taylor coefficients f^^^KO),f^^^-^\0) 
for each new value of j (the degree of the singularity). This is the principal obstacle to 
overcome. 

We overcome it in [|6] as follows: The expression ([9]) for the Balian-Bloch invariants of 
7, 7^, . . . consists of two types of terms, in terms of their dependence on the iterate r. They 
have a common factor of 2rL{hliy~''^Ar{0), and after factoring it out we obtain one term 

2 — 2 cos a/2 

which depends on the iterate r through the coefficient {hHY, and one 

Ki-,V-) (E(^2^)'/^'Ho)/(^^-^no)j 

which depends on r through the cubic sums Ylq=i{h2r)^ of inverse Hessian matrix elements 
/igr- In order to 'decouple' the even and odd derivatives, it suffices to show that the functions 
{hll)"^ and X]qLi(^2r)^ ^^6, at least for 'most' Floquet angles a, linearly independent as 
functions of r G Z, i.e. that (hH)^'^ YlgLii^lr)^ is a non-constant function of r. It is 
convenient to use the parameter a = — 2cos| and we write the dependence as h^j-^a). 
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We therefore define tlie 'bad' set of Floquet angles by 

(10) B = {a : the sequence {(^2r('^))~^ X^(^2r('^))^' r = 1, 2, 3, . . . } is constant in r}. 

9=1 

Using facts about the finite Fourier transform and circulant matrices, we compute that B = 
{0, 1, ±2}. Since the proof is computational, we also present a simple conceputal argument 
(cf. Proposition 16. 71) that B is finite, although the proof only gives the poor estimate 3^° 
on its number of elements. For Floquet angles outside of B, we can determine all Taylor 
coefficients f^^ (0) from the wave invariants and hence the analytic domain. 

We use a similar strategy in the dihedral D„-case in ^|71 Due to the extra symmetries, the 
inverse results in the dihedral case require much less information about the wave invariants 
than in the one symmetry case. 

1.3. Related results, (i) We have already mentioned the prior result that analytic drum- 
heads with up/down and left /right symmetries are spectrally determined in that class \Z1\ 
IZ2j . Previously, it was proved by Colin de Verdiere |CV] that such domains are spectrally 
rigid. To our knowledge, the only other prior result giving a 'large' class of spectrally do- 
mains is that of Marvizi- Melrose [MMlj . in which members of a spectrally determined two- 
parameter family of convex plane domains are determined among generic convex domains 
by their spectra. 

(11) In |Z4j . we extend the inverse result to the exterior problem of determining a Z2- 
symmetric configuration of analytic obstacles from its scattering phase (or resonance poles). 
Our result may be stated as follows: Let = _ [OUt^^lO} where O is a convex analytic 
obstacle, where x G O and where t^^l is the mirror refiection across the orthogonal line 
segment of length L from x. Thus, {O U Tx^l{0)} is a Z2-symmetric obstacle consisting of 
two components. Let denote the Dirichlet Laplacian on VL. We have: 

Theorem 1.5. [Z4J With the same genericity assumptions as in Theorem \l.l\ the resonance 
poles of determine O within the class of Z2 symmetric analytic obstacles. 

1.4. Future directions. An obvious future direction is to study the wave invariants with- 
out any symmetry assumptions. As will become clear from the calculations in this article 
(cf. Theorems 14.21 and 13. II) . symmetries make 'lower order derivative data' in wave invariants 
redundant and allow one to concentrate on terms in a given wave invariant with maximal 
numbers of derivatives. Lacking symmetries, the lower order derivative data is no longer 
redundant and one has to navigate a complicated jungle of terms to determine which com- 
binations are spectral invariants. It is plausible that one cannot work with just one orbit 
but must combine information from two bouncing ball orbits (they always exist in a convex 
plane domain). The main problem is then to extract from the wave invariants of the iterates 
of each bouncing ball orbit sufficient Taylor series data at the endpoints to determine the 
domain. To do this, it seems necessary to analyze how Feynman amplitudes of labelled 
diagrams behave as a function of the iterate r of the orbits. The graphs themselves do not 
depend on r, so the dependence comes from the labelling. 
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2. Billiards and the length functional 

We begin by establishing notation on plane billiards and length functions. After recalling 
basic notions, we calculate the Hessian of the length functional at iterates of a critical 
bouncing ball orbit in Cartesian coordinates adpated to the orbit. 

We denote by a simply connected analytic plane domain with boundary dQ of length 
27r. The billiard flow $* of Q is the broken geodesic of the Euclidean metric on Q. That is, 
for (x,^) G T*Q°, the trajectory ^) follows the Euclidean straight line in the interior 

Q° of fl and reflects from the boundary by Snell's law of equal angles. By the billiard map 
f3 oi Q we mean the map on B*dQ induced by $*: we add a multiple of the inward unit 
normal Ug to (g, rj) G B*{dQ) to obtain an inward pointing unit vector v at q. We then follow 
the bilhard trajectory t>) until it hits the boundary, and then define P{q,ri) to be its 
tangential projection. We refer to [PSl IKTl [Z5] for details and discussions of the billiard 
fiow on domains in M^. 

It is natural at first to parametrize dQ by arclength, 

(11) g : T ^ 9fi C M^ 

starting at some point go ^ 9Q. Here, T = M\27rZ denotes the unit circle. By an m-link 
periodic reflecting ray of Q we mean a periodic billiard trajectory 7 which intersects dQ 
transversally at m points q{^pi), ■ ■ ■ ,q{ipm), and refiects off dQ at each point according to 
Snell's law 

^ ^ ki^Hi) - qM\ ■ '^''^ \q{^j) - q{v,-i)\ ' 

Here, t'g(^) is the inward unit normal to dQ at q{(p). We refer to the segments q{(Pj) — q{(pj-i) 
as the links of the trajectory. We denote the acute angle between the link q{{pj+i) — q{ipj) and 
the inward unit normal i^g(<^^) by Z(g((y9j+i) — q{(Pj), ^q{<pj)) and that between q{(pj) — q{ipj^i) 
and the inward unit normal at q{fj) by Z{q{(pj) — g(y9j_i), t'g(^j)), i.e. we put 



For notational simplicity we often do not distinguish between a billiard trajectory in S'*f2 
and its projection to f2. 

We define the length functional on T*^ by: 

(14) l((/?i,...,¥?m) = \q{}p\) -qi}P2)\ + ■■■ + \q{}PM-x) -g(¥^Af)| + \q{}PM) -qi}P\)\- 

We often use cyclic index notation where g(v9M+i) = qis^x)- It is clear that L is a smooth 
function away from the 'large diagonals' Aj^+i := = v^j+i}? where it has singularities. 
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We have: 



(15) - = smZ(g((^j+i) - q{<^j), Ug(^^)) 

Hence, the condition that -^L = is the same as f|T2|) for the 2-hnk defined by the triplet 

Let 7 denote a periodic refiecting ray of Q. The hnear Poincare map of 7 is the 
derivative at 7(0) of the first return map to a transversal to at 7(0). By a non-degenerate 
periodic refiecting ray 7 we mean one whose linear Poincare map has no eigenvalue equal 
to one (cf. |PS[ IKT] ). The following relates P^ and the Hessian of the length functional in 
angular coordinates: 



Proposition 2.1. ( [KT] (Theorem 3)) Let denote the Hessian of L in angular coordi- 

det(/ -P^) = - det{-H^) ■ (61 ■ ■ ■ b^)-\ 



nates cpj at a critical point 7, and let bj = ^ ^'^^q^^q^ ^ Then 



This identity may be proved by expressing both sides in terms of bases of horizontal and 
vertical Jacobi fields. 

2.1. Cartesian coordinates around bouncing ball orbits. We now specialize to the 
case where 7 is a bouncing ball orbits (i.e. 2-link periodic refiecting rays). As in the 
Introduction, we orient fl so that the bouncing ball orbit is along the y-axis with endpoints 
A = (0, -j), B = (0, — -j) and parametrize d^l near Ahj y = f+{x) and near B hj y = f-{x). 
We do not assume the domain is up-down symmetric. 

We denote by Ra, resp. Rb, the radius of curvature of Q at the endpoints A, B. When 7 
is elliptic, the eigenvalues of P-y are of the form {e^*"} {a G M) while in the hyperbolic case 
they are of the form {e^"} (a G M). They are given by the same formulae in both elliptic 
and hyperbolic cases: 



(16) 



cos(a/2) = (elliptic case) 



cosh (a/2) = i^), (hyperbolic case) 



We define the length functionals in Cartesian coordinates for the two possible orientations 
of the rth iterate of a bouncing ball orbit by 



2r 



(17) C±{Xi, . . .,X2r) = X] V ~ ^j)^ + ifw±(3+l)iXj+i) - fw±(j)iXj)y 



Here, w± : Ij2r where w+{j) (resp. w_(j)) alternates sign starting with w+(l) = + 

(resp. if_(l) = — ). Also, we use cylic index notation where X2r+i = Xi. 
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We have: 



dxi 



(18) 



(xj-a;^-+l)+(/„^(j)(xj)-/„^(^+l)(a;j+i))/^^j^j(xj) 
^(xj-Xj+i)2+(/„^(j)(a;j)-/„^(j4.i)(xj+i))2 

(^■j-l -^j)+{/»± 0-1) 0) (j) 

A/(a^j-a^j-i)^ + (/u.±(j)(a^j)-/u.±0-i)(2^j-i))^ 



We will need formulae for the entries of the Hessian of at its critical point (xi, . . . , X2r) 
in Cartesian coordinates corresponding to the rth repetition of a bouncing ball orbit. 



Proposition 2.2. Put 



L 



a = -2(1 + Lfim = -2(1 -j^),b= -2(1 - Lfim = -2(1 - — ). 

Then the Hessian of at x = in Cartesian graph coordinates has the form H2 
/a 2\ 

-J- I for r = 1 and for r > 2, 

\2 b 



2r 



-1 

T 



1 







1 a 1 

16 1 

1 ... b 



Proof. A routine calculation gives 



m 



2(i + ^4j)/:,o)(0)), 



-(0) = X 

for r > 2. In the case of r = 1, the length functional is 2||(xi, /+(a;i)) — (x2, /-(2;2))||- Note 
that for r > 2, there are two terms of contributing to each diagonal matrix element and 
one to each off-diagonal element, accounting for the additional factor of 2 in the diagonal 
terms. Also note that fw+(j){0) - fw+(j+i){0) = w+{j)L and that /^(O) = f'l{0) = 

□ 

We remark that the Hessian in Cartesian coordinates in Proposition l2. 21 differs from that in 
angular coordinates in |KTj in that the off-diagonal entries differ in sign. This is because the 
graph parametrization gives the opposite orientation to the tangent TaOQ than the angular 
parametrization and the same orientation at TbOQ. The angular Hessian is related to 
the Cartesian Hessian by H^.^ = JH2rJ^ where J = diag{l, — 1,1,— 1,.. .,1,-1) is the 
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change of basis matrix. Clearly, the determinants of the two Hessians agree. Since bj = 
we obtain from Proposition 12.11 the following: 

Corollary 2.3. As above, let denote the Hessian of in Cartesian coordinates at 
the rth iterate Y' of a bouncing ball orbit 7 of length 2L. Then 

det(J -Py) = -L^" det{H2r). 

The determinant det is a polynomial in cos ^ (elliptic case), resp. cosh ^ (hyperbolic 
case) of degree 2r. In the following we restrict to the elliptic case. 

Proposition 2.4. We have 

det H2r = -L-"^'' (2-2 cos ra) . 

Proof. Let A^, be the eigenvalues of P^r, so that det{I — P^) = 2 — (A.,. + A^^). Now, if 
the eigenvalues of P^ are {e^*°} (in the elliptic case) then those of P^ are {e^*^"}, hence 
det{I — P^r) = 2 — 2cosra. Similarly for the hyperbolic case. The formulae then follows 
form Corollary 12. 3[ 

□ 

We now consider the inverse Hessian ?i+ = H2r , which will be important in the calculation 
of wave invariants. We denote its matrix elements by hF^ . We also denote by Ti- the matrix 
in which the roles of a, b are interchanged; it is the inverse Hessian of 

Proposition 2.5. The diagonal matrix elements are constant when the parity of p is 
fixed, and we have: 

p odd =^ = h^, p even =^ = h"^ 
p odd =^ h^]' = h^_^, p even =^ h^J' = h"^, 

Proof. Indeed, let us introduce the cyclic shift operator on M?"^ given by Pcj = e^+i, where 
{cj} is the standard basis, and where Pe2r = ei. It is then easy to check that PH+P'^ = H-, 
hence that PH^^p-^ = nz\ Since P is unitary, this says 

h^.' = {n-Je„e,) = {Pn^'p-'e„e,) = {H^'P^'e,, p-\) = h^-'^'^-'. 

It follows that the matrix T-C± is invariant under even powers of the shift operator, which 
shifts the indices j ^ j + 2k {k = 1, . . . ,r). Hence, diagonal matrix elements of like parity 
are equal. □ 

3. Resolvent trace invariants 

We now formulate the key results (Theorems 14. 2114. 2p expressing localized wave traces as 
oscillatory integrals over the boundary with special phases and amplitudes. We then tie 
these statements together with the statements in Theorem 1.1 (v) of |Z5] . 

First, we state a general result, largely contained in |Z4l IZ5j which expresses the localized 
resolvent trace as a finite sum of special oscillatory integrals. For simplicity we only state it 
for the rth iterate of a bouncing ball orbit. 
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Theorem 3.1. Suppose that rLj is the only length in the support of p. Then for each order 
in the trace expansion of Corollary (3^, we have 



± M:2r<M<R+2r u:\u\<R,M-W\='ir 

where a runs over all maps a : {1, ...,M} — > {0,1}, and where l^M^ik^ are oscillatory 
integrals of the form 

(19) 



X 0''^4"p{k, Xi, X2, . . . , X2r)dXi ■ ■ ■ dx 



2r- 



Here, C^,^ is given in ( [77| ) and a'^^ are certain semi-classical amplitudes (cf. ^^)- The 
asymptotics are negligible unless M —\a\= 2r and then the order of I'^^^{k) equals — |cr|. 

It follows that only a finite number of terms I'^^j"^{k) contribute to each order in k in the 
expansion in Corollary | 



Corollary 3.2. We have: 

where B^.j are the Balian-Bloch invariants of the union of the periodic orbits 7, and 'DB,-y{k + 
it) is the symplectic pre- factor of 

3.1. Proof of Theorem l3.1i As mentioned above, most of the proof is contained in [Z4l[Z5]. 

For the sake of completeness, we sketch the key elements of the proof. 

We follow the path originated by Balian-Bloch and followed in many physics articles (see 
e.g. [BBTI [BB2| IAG] ). It starts from the exact formula (of Fredholm-Neumann) , 

(20) R%{k + ir) = Ro{k + it) - 2 Ve{k + it){I + N{k + ir))-^rnSt{k + ir) 

for the resolvent with given boundary conditions. Here, T>^{k + it) (resp. Sl{k + ir)) is 
the double (resp. single) layer potential, S^'^ik + ir) is the transpose, and N{k + it) is 
the boundary integral operator on L'^{dQ) induced by Vi{k + ir). Also, Ro{k + ir) is the 
free resolvent on M^, and is the restriction to the boundary. The Schwartz kernel of the 
boundary integral operator is given by plus (in the Dirichlet case) or minus (in the Neumann 
case) 

(21) N{k + iT)f{q) = 2 [ -^Goik + IT, q, q')f{q')ds{q'), 

J an 

where Go{X,x,y) is the free Green's function (resolvent kernel) on M^, where ds{q) is the 
arc- length measure on dQ, where z/ is the interior unit normal to Q, and where = ■ V. 
The free Green's kernel has an exact formula in terms of Hankel functions fl3ip . which gives 
a WKB approximation to N{k -\- ir) away from the diagonal. Its phase is the boundary 
distance function d^{q, q'), indicating that N{k + ir) is the quantization of the billiard map. 

But as discussed extensively in \Z5\ [Z4t IHZj . N{k + ir) is not a classical Fourier integral 
operator, but is rather a non-standard kind of hybrid Fourier integral operator. Near the 
diagonal, it is a homogeneous pseudo-differential operator of order —1 (in dimension two it 
is actually of order —2 as proved in [Z5J, Proposition 4.1), while away from the diagonal 
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it is a semi-classical Fourier integral operator of order which quantizes the billiard map. 
To separate out these two Lagrangian submanifolds (which intersect along tangent vectors 
to the boundary), we introduce a cutoff xi^^'^ll ~ q'\) fo fhe diagonal, where 6 > 1/2 and 
where x ^ C*^(]R) is a cutoff to a neighborhood of 0. We then put 

(22) N{k + iT) = Noik + ir) + Ni{k + iT), with 



(23) 



No{k + iT, q, q') = x{k^'^\q " q'\) N{k + ir, q, q'), 



N,{k + IT, q, q') = (1 - x{k^-'\q " q'\)) N{k + it, q, q') 



As proved in \Lb\ [Z4l IHZj . Ni{{k + it), q, q')) is a semiclassical Fourier integral operator of 
order with phase equal to the boundary distance function dg^lq, q'). The diagonal part A^q 
is of order —1 (in fact, of order —2 jZ5j ) and therefore plays a secondary role. 

We now relate the expansion of the regularized resolvent trace to that for logdet A^(A; + 
it). This relation has already been proved in |EP[ [U| IZ4j in somewhat different ways. 

The clearest proof is to combine the interior boundary problem with a complementary 
exterior boundary problem A^^ Since we are only dealing here with Dirichlet or Neumann 
boundary conditions, we do not define the term 'complementary' but only use the term to 
indicate the special cases B = D,B' = N or B = N,B' = D. We therefore introduce the 
exterior Green's kernel G%{k + iT,x,y) G V{Q'^ x Q'^) with boundary condition B, namely 
the kernel of the exterior resolvent and is the unique solution of the boundary problem: 



(24) 



-(Ag: + {k + iTf)G%{k + IT, X, y) = 5y{x), {x, y e Q') 
B'G%{k + iT,x,y) = Q, xedVl" 

—^—Q-^ i[k+iT)G'^,[k+iT,x,y) = o[-), as r oo. 



We now combine the interior and exterior operators with complementary boundary con- 
ditions B, B' into the direct sum R%{k + it) © R^i{k + it). For simplicity, we only consider 
B = D,B' = N. For p G C;^(M+), we put 

(25) R%{k + iT) © R^B, {k + iT) = p{k - fi){fi + iT) [R% {fi + ir) © Rf, {fi + iT)] dp. 

The purpose of combining the interior /exterior resolvents is revealed in the following propo- 
sition, which equates the trace of the direct sum resolvent to the Fredholm determinant of 
the boundary integral operator. It is proved in jZ4j and closely related statements are proved 
in [EP\ \C\ . The operator A^ is defined in 0211) in the Dirichlet case. In general it depends on 
the boundary conditions B, B' . We follow the notation of [Tj except that we multiply the 
A^ of [T] by i to simplify some notation. 

Proposition 3.3. For any r > 0, the operator {I + N{k + iT)) has a well-defined Frehdolm 
determinant det(/ + A^(A + iT)), and we have: 

Tr^2 [R^^{k + it) © R^;,ik + it) - R^p{k + it)] 

= j^p{k-\)j^\ogdei{I + N{\ + iT))d\. 
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Further, for r > 0, logdet(/ + N{k + ir)) is differentiable in k, (J + N{k + ir)) ^N'{k + ir) 
is of trace class and we have: 

4r^ogdet(I + N(k + iT)) = TraM + Nik + iT))-^N'{k + it). 
dk 

This proposition reduces wave trace expansions to the boundary. Indeed, the direct sum 
resolvent is related to the direct sum wave groups as in ([7]): 

/"OO 

(26) R%{k + it) © R%,{k + it)= / p(t)e*('=+^^)* © dt. 

Jo 

The trace of the direct sum wave group (t) © E^i (t) has a singularity expansion as in (jlj) 
which sums over interior and exterior periodic orbits. As in ([8]), it may be restated in terms 
of the direct sum resolvent: Let 7 be a non-degenerate interior billiard trajectory whose 
length is isolated and of multiplicity one in Lsp{Q). Let p G C^{L^ — e,L^ + e), equal 
to one on {L^ — e/2, + e/2) and with no other lengths in its support. Then the interior 
trace TrR^p{k + it) and the exterior trace Tr[R^,p{k + it) — Rop{k + ir)] admit complete 
asymptotic expansions of the form 



(27) 



Tr[i?g:^(fc + it) - Ro,{k + it)] ~ VB,.y{k + it) J^JL^ 5^,, k- 
TrR%{k + it) ~ VB,-,{k + it) Y.7=^ ^"^ 



whose coefficients B^.j are the Balian-Bloch resolvent trace invariants of periodic (internal, 
resp. external) billiard orbits. We can therefore sum the two expansions to produce one for 
the direct sum. The coefficients depend on the choice of boundary condition but we do not 
indicate this in the notation. 
Combining the results, we get: 

Corollary 3.4. Suppose that is the only length in the support of p. Then, 
f^p{k - A)^ logdet(/ + iV(A + tT))dX 
= J^p{k-\)Tran{I + Ni\ + iT))-^N'{X + iT)dX , 

~ 'DB,yik + zr) Y.f=o k~^ 

where as above B^.j are the Balian-Bloch invariants of the union of the periodic orbits 7 of 
length of the interior and exterior problems in ( [F7| j. 

In proving the remainder estimate and the expansion in Proposition 13. 6^ we further mi- 
crolocalize the result to the (interior) orbit 7. This will select out the wave invariants of the 
desired interior orbit 7. A periodic orbit of the billiard flow corresponds to a periodic point 
of the billiard map (3. To microlocalize to this periodic orbit we introduce a semiclassical 
pseudo differential cutoff operator Xo(v^, k~^D^). In the case of a bouncing ball orbit, it has 
complete symbol xiViV) supported in K := {(v^,??) : Iv^l, \ri\ < e}. 
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Proposition 3.5. Suppose that is a bouncing ball orbit, whose length Lj is the only length 
in the support of p. Let xo be a cutoff operator to the endpoints of'-f. Then, 

Trp * (/ + N{k + ir))-^ o -^N{k + ir) 

~ Trp * (/ + N{k + zr))-i o ^N{k + ir) o xo{k). 

We will use the formula in Corollary 13.41 as modified in Proposition 13.51 to calculate the 
B^-j modulo remainders which are inessential for the inverse spectral problem. To do so, we 
now express the left hand side (for each order of singularity k~^) as a finite sum of oscillatory 
integrals /^'^ (see (|T9|) ) plus a remainder which is of lower order than k~K 

To define the oscillatory integrals I^}'^, we first expand (/ + A^(A + «t))~^ in a finite 
geometric series plus remainder, 

Mo 

(28) (/+iV(A + ?r))-i = 5^(-l)''' N{\ + irf' + {-l f'^^+^ iV(A + zr)*^°+^(/+iV(A + zr))-\ 

M=0 

and prove that, in calculating a given order of Balian-Bloch invariant B^j, we may neglect 
a sufficiently high remainder. 

Proposition 3.6. For each order k~'^ in the trace expansion of Corollary ^3.^^ there exists 
Mq[J) such that 

(0 E;to(-l)*'^^ k P(k - A) N{X + zT^N'iX + zT)dX 

= VB,y{k + it) + 0{k-J-^), 

{ii) Tr J^p{k - X)N{X + iTy"°+\I+N{X + ir)y^N'{X + iT)dX = 0{k-^-^). 

The same holds after composition with Xo(^)- 

The proof of this Proposition is one of the principal results in |Z5l IZ4j . In |Z5] the result 
is stated in Theorem 1.1 (iii), while the remainder trace is estimated in §8. The version 
stated in Proposition 13.61 is proved in §5 of [Z4j. It is simpler than Theorem 1.1 (iii) of |Z5j 
because the interior integral analyzed in §7 of that paper is eliminated in the reduction to 
the boundary. 

It simplifies the formula somewhat to integrate the derivative by parts onto p, since it 
eliminates the derivative in the special factor A^'(A + ir). 

Corollary 3.7. For each order k~'' in the trace expansion of Corollary \3.4\) there exists 
Mo{J) such that 

EMlo^-i^TrJ^p'{k-X)N{X + ZTr^'dX 

= VB,,{k + ir) E/^o 5^,, k-^ + 0{k-J-'), 

{ii) Tr J^p{k - X)N{X + iTy"''+\I+N{X + ir)y^N'{X + iT)dX = 0{k-^-^). 
The same holds after composition with Xo{k)- 
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The next step is to prove that the terms in Proposition I3.6( i) may be expressed as oscil- 
latory integrals (see (fT9l) ). This is not obvious, as mentioned above, since the N operator is 
not a Fourier integral kernel. As indicated in (!22|) - (!23|) . we handle this problem by breaking 
up iV as a sum N = Nq + Ni of two terms, where Nq has the singularity on the diagonal of a 
pseudodifferntial operator of order —2 (cf. [Z5], Proposition 4.1), and where A^^i is manifestly 
an oscillatory integral operator of order with phase \q{<f) — As mentioned above, 

and as discussed in detail in [Z4l IHZj . the phase is a generating function of the billiard map, 
so the Ni term is a quantization of f3. 

We thus write, 

(29) {No + N,f'= Yl Ar,(i)°Ar,(2)0---oAr.(M). 

o-:{l,...,A/}^{0,l} 

In |Z5] §6, we regularized the terms by proving a composition law for products NqoNi, NioNq. 
The main technical point is that the amplitudes of A^^o, A^i belong to the symbol class >S'^(T) 
where T is the unit circle parameterizing dQ, consisting of symbols a{k, ip) which satisfy: 

(30) \{k-'D^ra{k,y^)\ < C,|A;r^l"l, (|A;| > 1). 
This follows from the classical formula (see e.g. |Z5] §4; |AGj . (2.2)) 

N{k + ^T,q{^^),q{^2)) = -^k + iT)H['\{k + irM^,) - q{^,)\) 

(31) 

X cos Z(g(y?2) - qM, z^g(^2)) 

for A^ in terms of Hankel functions and from the asymptotics of Hankel function h[^\ We 
recall that the Hankel function of index u has the integral representations ([TJ, Chapter 3, 
§6) 

(32) Hi'\z) = (^)^/^ "';;r;;r - i-.r-'^'ds, 

from which it follows that h[^^ admits an asymptotic expansion as its argument tends to 
infinity of the form 

oo 

(33) H[^\t) ~ e'*-^r^ Yl ^i^"'' ^ 

j=0 



where cq = A/2/7r. Moreover, the expansion can be differentiated term by term. We set: 



2 -oo 



(34) '-^"'^'^ - i'"^*- 

SO that 

(35) H[^\t) r^e''-''^t--2a^{t) 

We observe that ai is a complex valued semi-classical symbol of order of 2; G IR+ in the 
sense that (cf. ( l30l) ) 

(1 - xik'-'z))ai{ik + tT)z) e 5°(M,). 
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We then have 

(36) {k + iT)H[^\{k + iT)z) = {t^)^e'^''^^^'ai{{k+tT)z), 

hence 



N,{k + IT, q{ip,), g(^2)) = (1 - Xik'-'ivi - ^2))) 
i \,^^')tl^,)\ y^a,{k + tT,q{^,),qM^^^ 

with 

(37) a,{k + IT, g((^i), q{^2)) := a,{{k + ir)|g((^i) - q{^2)\) co^^i,2 G ^^(T^), 

where ^^1,2 = ^q{V2) - q{Vi), ^q{ip2))- 

The main conclusion is that NqNi and NiNq are semiclassical Fourier integral operators 
with the same phase as A^^i, but with an amplitude of one lower degree in k. This allowed us 
to remove all of the factors of A'o from each of these terms except for the term Nq^ . Each 
remaining term except for Nq^ is a Fourier integral operator on for some m < M, with 
phase given by the length functional ([T^ and with amplitude in the symbol class S^{T"^) 
for some p, consisting of symbols a{k, ipi, . . . , ip^) which satisfy the analogue of fl5U]l : 

(38) \{k-^D^Ya{k,^)\<CM^~'\''\ (|fc|>l). 

Because each removal of A^o drops the order by one, the term Nf'^ is of the highest order 
in the sum. A later estimate on traces shows that Nq '^ does not contribute to the trace 
asymptotics (see |Z5j, §9.0.7). 

We summarize the result as follows. Let us rewrite the terms of fl29l) as 

(39) N„ := Ar,(i) o Ar,(2) o • • • o N^^m) 
and set 

(40) |cr| = ^a^'^(0) = the number of A^^o factors occurring in No-. 

In |Z5j . Propositions 6.1, we show that the regularized compositions are semiclassical Fourier 



integral kernels. 
Proposition 3.8. We have: 

(A) Suppose that is not of the form N^'^ . Then for any integer R > 0, o xoik + ^t) 
may he expressed as the sum 

No = Fo{k,ipi,ip2) + Kn, 

where is a semiclassical Fourier integral kernel of order — associated to of the 

form 

(41) Fo{k, if,, ^2) = e^('=+-)l'?(^i)-.(^2)|^^(^^ 

where Aa-{k,ipi,ip2) is a semi-classical amplitude, and where the remainder is a bounded 
smooth kernel which is uniformly of order k~^. 

(B) A"^ o Xo ~ Nqm o X05 where Nqm is a semiclassical pseudodifferential operator of order 
—M. (For the notation xo see Proposition \3.5^) .) 
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As a corollary of Proposition 13.81 we obtain the following preliminary form for the trace 
as a sum of oscillatory integrals. It is a simplification of |Z5] . Lemma 9.2 in that we do not 
need any interior integrals. 

Corollary 3.9. Tr p' * o is an oscillatory integral of the form 

rcr (U\ _ u(M-\a\+3)/2 f f f ik[{l-ti)t+tiC^{q{ipi),...,q{ipM^\„\))] f,~rlogkC^{q{ipi),...,q{ipM_^„^)) 



where xili^i) ~~ 1(^2), V^i) is the value at the vector {q{(pi), q{(pi) — q{^2)) of a cutoff x to 
a microlocal neighborhood in B*dQ of the direction of the bouncing ball orbit, where 

/:^(g(V9i),...,g((^A/-|a|)) = \qM - qMI + ■■■ + \q{VM-\a\) 

and where Al,j{k, (pi, . . . , v^Af-io-i) £ 5*^'°^'. 

3.1.1. Completion of proof of Theorem \3.1[ We now complete the proof of Theorem 13. II To 
obtain our final form for the oscillatory integrals, we make some further simplifications. For 
simplicity of exposition, and because it is our main application, we specialize to a bouncing 
ball orbit. In view of Propositions 13.31 and 13. 6[ it suffices to prove: 

Proposition 3.10. Suppose that rL^ is the only length in the support of p. Then for each 
order k~^ in the trace expansion of Corollary ([J.^p, we have 



[ p{k~X)-^\ogdet{I + N{X + ^T))dXr^Yl E E ^m7(^) 



± M:2r<M<R+2r a:\cr\<R,M-\cr\=2r 

where the oscillatory integrals I'^fp^{k) are as in Theorem \3.1[ 

Proof. The first observation is that the regularized integral I'^.j p{k + iT) of Corollary 13.91 has 
no critical points unless M — \a\ = 2r (where rL^ is the unique length in the support of p). 
We will refer to these oscillatory integrals as contributing. Since each has two pieces, each 
contributing integral can be written as a sum of 2^'' terms I'^p{k + zr), corresponding to a 
choice of an element w of 

{±Y^ := {w : Z,r - {±}}. 
The length functional in Cartesian coordinates for a given assignment w of signs is given by 

2r 



(42) C^ixu . . .,X2r) = E V (^i+i-a:i)H(/^(jH-i)(a;j4_i)-/^(j)(a;j))2. 

i=i 

Here, X2r+i = Xi. 

We further observe that Ip^p{k + ir) has no critical points unless w{j) alternates between 
+ and — as j increases. Otherwise, I^^^pik + ir) is negligible as A; — > 00. Thus, only two 
w count asymptotically, which we denote by w±. The corresponding length functionals are 
given in (ITSl) and their Hessians are given in Proposition 12. 2[ 

In these remaining oscillatory integrals, we then eliminate the (t, p) variables in the integral 
displayed in Corollary 13.91 by stationary phase. The Hessian in these variables is easily seen 
to be non-degenerate, and the Hessian operator equals — ^q^q^ ■ The amplitude depends on 

t only in the factor p'{t). Since p'(t) = tp{t) and since p is assumed to be constant in some 
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interval {rL^ — e,rL^ + e), tp{t) is locally linear and therefore only the zeroth order and 
(— l)st order terms 



ik dk 

in the stationary phase expansion are non-zero. In the second term, the k in the denominator 
comes from the Hessian operator and the k in the numerator comes from the /i- derivative 
of the amplitude. After replacing the dtdfi integral by this stationary phase expansion, we 
arrive at the final form of the oscillatory integrals (fT9|) given in the Theorem, with amplitude 

(43) air{k,x) = C^,Al,{k,x) + ] ^'^\f''\ k,x). 

□ 



4. Principal term of the Balian-Bloch trace 

In this section, we state and begin the proof of a key result for the proof of Theorems 
11.11 and II. 4[ It singles out a single oscillatory integral (the principal term) from Theorem 
13.11 which generates all terms of the wave trace (or Balian-Bloch) expansion which contain 
maximal number of derivatives of the boundary defining function per power of k (i.e. order 
of wave invariant). As mentioned in the introduction, the other terms will turn out to be 
redundant for domains in our sjnnmetry classes. 

To clarify this notion of generating all the highest derivative terms, we define it formally. 
Below, denotes the s-jet. 

Definition 4.1. Let 7 be an m-link periodic reflecting ray, and let p G C^(M) be a cut off 
satisfying supp p{~\Lsp{VL) = {rL^} for some fixed r E N. Given an oscillatory integral I{k), 
we write 

TrlnRspik + tr) = I{k) mod C(^ k-^{J^^-^K)) 

j 

^f 

Trl^R%{k + iT)-I{k) 

has a complete asymptotic expansion of the form and if the coefficient of k~^ depends 
on < 2j — 2 derivatives of the curvature k at the reflection points. 

For the sake of clarity, we state the next result only in the simplest case of a bouncing ball 
orbit. The statement is similar for any non-degenerate m-link periodic reflecting ray. The 
description of the properties of phase and amplitude are repeated from |Z4j for the sake of 
self-completeness. For terminology concerning billiard trajectories, we refer to ^ 

Theorem 4.2. Let be a primitive non-degenerate 2-link periodic reflecting ray, whose 
reflection points are points of non-zero curvature of dQ, and let p G C^(R) be a cut off 
satisfying supp p fl Lsp{Q) = {rL^} for some fixed r G N. Orient Q so that 7 is the vertical 
segment {x = 0} (1 fl, and so that dQ is a union of two graphs over [— e, e]. Then in the 
sense of Definition \4-l\ we have 
(44) 

TrlnR%{k + tr) = E± I[-e,e]^r e*('=+-)^±(-i---2.)p(£^(a;^^ . . . ^ X2rWl,{k, x,, X2, . . . , X2r)dXi 
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where the phase . . . ,X2r) is given in ( (77| j, and where the amplitude is given by: 



1 d 

'^±,ri^^ . . . , X2r) = ^w±A±^ri^i . . . , X2r) + ' ' ' ' *''2r); 



where 



(45) 



xpr /7 \ _ T-r2r , ai ((fc+jr) ^ {xp-Xp^i)H{f^^ (^j {a;p)-/„^ (j^i) (xp+i )) 



X 



{Xp-Xp+l)f!^^^^^ {Xp)-{f^^(^p) {Xp)-f^j^(^p+l',{Xp+l))^ 



y/{Xp-Xp+l)2 + {f^^(^p^{Xp)-f^^^p+l^{Xp+l))2 

where ai is the Hankel amplitude in ( f^) . Here, as above, X2r+i = xi. 



Theorem 14.21 is a crucial ingredient in the proof of Theorem ll.li It gives exphcit formulae 
for the phase and amplitude of the principal oscillatory integrals that determine the highest 
order jet of fl in each wave invariant. The notation A^"^, a^^ refers to the amplitude of the 
principal terms of the 2rth integral; these amplitudes contain terms of all orders in k and 
principal here does not refer to the principal symbol, i.e. the leading order term in the 
semi-classical expansion. The calculation of the highest derivative terms of the Balian-Bloch 
wave invariants uses only some key properties of the phase and principal amplitude which 
may be derived directly from the formulae in Theorem 14. 2[ They are detailed in §4.11 

The proof of theorem 14.21 requires two main steps: 

(1) Identification of two main terms in Theorem 13. 11 the principal terms, which generate 
the highest derivative data, and proof that the amplitude and phase have the stated 
form. 

(2) Proof that non-principal terms contribute only lower order derivative data. 

We now define the principal terms. In §4.H Lemma 14.51 we prove that their phases 
and amplitudes have the stated form. We further describe the properties of the phase and 
amplitude which will be used in the proof of Theorem 11.11 and tie the statement of Theorem 
14.21 together with the corresponding statement in [Z5] . The fact that non-principal terms do 
not contribute highest order derivative data to a given Balian-Bloch invariant requires the 
analysis of the stationary phase expansions in the next section and is given in §5.41 

Definition 4.3. Let 'j be a 2-link periodic orbit. The principal terms are the completely 
regular terms l2r'^^ coming from Nf^ ,i.e. with M = 2r and with ao(j) = 1 for all j . The two 
terms correspond to the two possible orientations w±{j), of the 2rth iterate of the bouncing 
ball orbit. 

In other words, the principal terms are simply those coming from the term 
(46) Trp* Nl'\k) o N[{k) o x{k) 



in the expansion 

We observe that in fact, the two principal terms are equal. This is not surprising, since a 
bouncing ball orbit is reciprocal. 

Proposition 4.4. We have: /2?,p'"+(A;) = /2?,p'""(A;). 
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Proof. We permute the variables xj according to the cychc permutation s of their indices: 

/ 1 2 ■■■ r-1 r \ 

\ 2 3 ■■• r 1 J 

in the integral in f|T9|) . Since = this takes £_ — > and a° ^ in fHSjl . 

Indeed, £± (resp. a^.) are sums (resp. products) of terms of the form F{xp—Xp+i, fw±{p){xp) — 
fw±{p+i){xp+i)). Cyclically shifting the index by one moves each term (resp. factor) to the 
next except that it does change the index w±{p). Hence, it changes the sum (resp. product) 
only by shifting to (and vice- versa). 

□ 

Henceforth, we often omit l2r,'p~i^) multiply I^r^^^ik) by 2. 

4.1. Key properties of the principal amplitude and phase. We first prove that the 
phase and amplitude of the principal oscillatory integrals have the form stated in Theorem 
14.21 and establish a few consequences. After that, we assemble all of the properties used in 
the proof of Theorem 11.11 In the following, we abbreviate £+ = C^,^ . We use the notation 
Dj-p = and multi-index notation for its powers. 

Lemma 4.5. The phase and principal amplitude of the principal oscillatory integrals 1^^'^^ 
have the following properties: 

[i) In its dependence on the boundary defining functions f±, the amplitude has the form 
ar{k,xJ±J'^). 

(a) As above, in its dependence on x 

A^_^^^{k, xi, . . . , X2r) = H2!l^Ap(xp, Xp+i) (2r + 1 = 1) 

{Hi) At the critical point, the principal amplitude has the asymptotics 

a^j.{k, 0) ~ {2rL)L~^Ar{0) + 0{k^^), where ^^(0) depends only on r and not on Q; 



(ma) "-^-^''"^^ ^detHes^z:, (^rL) A.(0) VB,,{k + zr)(l + 0{k-^)) (c/ilD; 

{iv)Va^j^^^{k,xi, . . . ,X2r)U=o = 0. 

{v) D(y-^)/:+U=o = 2^+(p)/i^^-^;)(0) mod R2riJ'^-'f+{0),J'^-'fA0)), 

(v.a) Di?')/:+U=o = 2wMfJlU0) modR2r{J'^-'fA0),J'^-^f40)), 
where = in general means equality modulo lower order derivatives of f . 
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Proof. The oscillatory integrals 12°'^^ have the form ( |T9i) with the phases C± (1421) . and by 
Proposition 14.41 it suffices to consider the + term. 

Formula (ii) for the amplitude follows from the general description of the amplitudes of 
all the oscillatory integrals in the proof of Theorem 13.11 (cf. fH3|) ). The factors A^^,. of 
the amplitudes of /2r,p"'* are given in fllSl) . 

The further properties of the phase and amplitude stated in Lemma 14.51 may be read off 
directly from the formula in ( H5l) . Statements (i)-(ii) are visible from the formula. At x = 0, 
the leading order term of the principal amplitude in k equals 2rL (from the factor C) times 
from the t~2 factor in the Hankel asymptotics times a coefficient ^^(0) which 

depends on r but not on Q and which is due to additional factors in the asymptotics of the 

free Green's function Gq: namely, a product of 2r factors of y^e^ from the principal term 

of the Hankel amplitude ai (loc. cit.), factors of ^ in the relation between the free Green's 
function Gq and the Hankel function ( pTl) . factors of 2 in the relation of N{k + it) and Go 
( I2TI) . We do not need to know ^,.(0) or other universal factors explicitly, since they multiply 
all terms in the expansion. Statement (iiia) gives the principal term in the stationary phase 
expansion at x = and relates the Hessian determinant and L"^' to the Poincare determinant 
as in Propositions 12.11 and 12.41 (see also [AGj . (3.17)). The second term is of order k~^, so 
will not contribute to the highest derivative term in a given wave invariant. 
From the fact that a; = is a critical point of f± and [xj — x,_i)^ we get 



(47) 



(a) ( ^ (xp - Xpi-i ) 2 + ( (p) (xp ) - {pH){x-M-iy)) \x=o = 

^ / {^■p-^'p+l)/^^(p){^p)-(/«J±(p)(^'p)-/™±(p+l)(^P+l))\ I Q 

^ I iy(xp-Xp+i)2+(/„^(j,)(a;p)-/„^(p+i)(xp+i))2 J 1^-0 



which implies 

(48) a+,rU=o = Va; Dka^_^,.\^=Q = 0. 

Statement (v) on the phase holds because 
(49) 

X {U+{p){xp) - /^+ (p±i) (xp±i ))/^^^^(pJ^(xp) 1^=0 mod i?2r(J^^"V±(0) ), 

D^x/^ C+\x=0 = X]±((^P ~ Xp±l)'^ + {fw+{p){Xp) — fw+{p±l){Xp±l)y)^^^'^ 

X {Uip)ixp) - Uip±i){xpii))f!;X)^Xp)U=o mod R2riJ''-'f+iO), J'^-'fAO))- 

We make the crucial observation that the ± terms are equal (and especially, do not cancel 
!), giving the factor of 2 in (v) since U+{p){0) - fw+ip±i){0) = w+{p)L. 

□ 



4.1.1. Further properties of the amplitude and phase. We continue the discussion of the 
amplitude by detailing the other special values of the phase and amplitude at the critical 
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point that are used in the ^ in the course of proving Theorem 11.11 Although the value of 
the discussion will only become clear in ^ it seems best to give the details at this point. 

(1) In the proof of Lemma [5.6( i). we use that 

(50) D^,-X'',rU=o = mod R2r{J''~'f±m, (Vp = 1, . . . , 2r). 

Indeed, by the explicit formula of (l45l) one can only obtain the higher derivative 
f'^~^{0) by applying all 2j — 2 derivatives on the term /^^(p)(a;p) in 

But then the accompanying factors of X2p — X2p+i vanish at the critical point. 

(2) In the proof of Lemma [5.6( ii). we use that 



(51) Df^-^^D,^C ^ mod R2r{J''-^ f±m. (Vp = 1, . . . , 2r, Vg ^ p). 

Indeed, in ( H9l) D^hJ'^^C is displayed as a product of two factors. Since q ^ p, the 
derivative D^^ must be applied to the factor 

{{Xp — Xp+i)^ + {fw+{p){Xp) — fw+{p+l){Xp+l))'^) ^^'^{fw+{p){Xp) — fw+{p+l){Xp+l))^ 

which vanishes at x = for any q. 
(3) In the same Lemma 15.61 we also use that the only non- vanishing third derivatives of 
£ at a; = are pure third derivatives in one variable DI..C Indeed, from (ITSl) . we 
see that only mixed derivatives using two consecutive indices (say, be 
non-zero. However, we have: 



(52) Dlp,^^,Cl=, = = D^^Dl C 



x=0- 



Since the identities are similar, we only consider the first, which is equivalent to 

'''''' - Xj+iy + ifw±{j)ixj) - fw±{j+i){xj+i)y 

We write the fraction as E^psEi+l^ and note that 

^x,^x,+i^U=o = ' U=o if F(0) = VG'(O) = 0. 

When F = ixj-Xj+i) + {f^^(j){xj)-f^^(j+i){xj+i))f^^f^.^{xj), we also have D^^D^^^^F\^=o 
0. 

(4) Further, we use that, for all p, Dl^C+{0) = 2w4.(p)/^'^(p-)(0). Indeed, as in the calcu- 
lation of the higher derivatives in Lemma 14.51 there are two terms, and each (in the 

notation above) has the form — q^qJ"^- obtain a non- zero term, the two derivatives 
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must fall on the factor /^^(p)(a^p), and thus we get 

D^JC+{0) = E±((%-a;p±i)^ + - /^+(p±i)(xp±i))2)-V2 

X {fw+{p){xp) - /«,+(p±i)(a^p±i))/lf_||(p)(a^p)U=o 

= 2^+(p)/S(P)(0)- 

Again, we observe that the Xp±i terms agree; therefore they add rather than cancel. 

4.2. Comparison with |Z5] . For the sake of completeness, we tie together the statement 
of Theorem 14.21 with the corresponding statement (v) of Theorem 1.1 of |Z5j and with [Z4] : 

Theorem 1.1 (v) of |Z5] .• Let 'y be a primitive non- degenerate m-link periodic reflecting ray 
of length L^, and let p G C^(M) be a cut off satisfying supp p fl Lsp{fl) = {rL^} for some 
fixed r G N. Then modulo an error term R2r{<J^'' '^ depending only on the (2j — 2)-jet 
of curvature n of dQ at the m reflection points aj ofj, the wave invariant B^j-i + B^-r j_i 
can be obtained by applying stationary phase to the oscillatory integral 

Tr p* A^f ^ o x{k) o Sl{k + irf o Vl{k + it). 



In Theorems 13.11 and 14.21 we have followed |Z4] in combining the interior and exterior 
problems. Taking the trace then eliminates the single and double layer potentials S^ resp. 
'D^ in Theorem l.l(v) of |Z5j . allowing for the reduction of the trace to the boundary in (H6|) . 

5. Feynman diagrams in inverse spectral theory 

In this section, we use the oscillatory integrals in Theorems l4.2l to obtain explicit formulae 
for the highest derivative terms of the wave trace invariants at a bouncing ball orbit in terms 
of the curvature function of the boundary. To our knowledge, these are the first explicit 
formulae. In the next section it will be proved that lower order derivative data is redundant 
for domains with our symmetries. 

For simplicity we restrict to bouncing ball orbits. There are similar results for general 
periodic reflecting rays (see Lemma [7.11 for the dihedral case). We first state the result for 
domains without symmetries, and then specialize to mirror symmetric domains in Corollary 
15.111 We use the graph parametrization rather than the curvature in the formulae. In the 
following, h^^ are the matrix elements of the inverse Hessian Hess{C+)~^ of the positively 
oriented length functional = C^^ of (ITS]) and in the principal terms. 

Theorem 5.1. Let Q be a smooth domain with a bouncing ball orbit 7 of length rL^. Then 
there exist polynomials P2,r,j{(,i, ■ ■ ■ , Gi+i! Vi^ ■ ■ ■ ^ V2j+i) which are homogeneous of degree —j 
under the dilation f — > A/, which are invariant under the substitutions C,j <^==4* —rjj and 
under f{x) f{—x) such that: 

. B,r, = P2,rAf^'\o), /l'^(o), ■ ■ ■ , f^y^'\oy, /f (0), /f (0), ■ ■ ■ , /f +'^(0)). 

• In the Balian-Bloch ( resolvent trace ) expansion of Corollary \3.4\ and in [2l\ ), the 
data f^''\o), f^'' ^\o) appear first in the k~^^^st order term, and then only in the 
expansion of the principal terms; 
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• This coefficient has the form 
B,r,_, ^ 4rLA(r){2(«; ..o)((/^ii,,,)Vf \0) - {h^.y f%m 

+R2r{J'^-'f+{0),J''-'f-m, 

where the remainder R2r{J^'^-'~'^f+{0), J''^^~'^f-.{0)) is a polynomial in the designated 
jet of f±. Here, = (—1)^'*'^ and as in the introduction, wg = \jiJt{g)\ '^'^^ 

combinatorial factors independent of Q andr. 

Where possible, we have simphfied the sums using Proposition [231 The top even derivative 
term is calculated in Lemma 15.51 and the top odd derivative is cacluated in Lemma I5.6[ 

The methods we use to make the calculations could be also used to evaluate the oscillatory- 
integrals in Theorem 13.11 and the wave invariants to all orders of derivatives. This could be 
useful in the inverse spectral problem for general domains without symmetry. However, we 
are content here to study the highest derivative terms and apply the results to domains with 
symmetry. 

We prove Theorem 15.11 by making a stationary phase analysis of the oscillatory integrals 
in Corollary 13. 1[ As mentioned in the introduction, our strategy involves a novel feature 
of the stationary phase expansion, namely to separate out the terms of the stationary each 
order in k which have the maximum number of derivatives of the boundary defining function 
or equivalently of its curvature. 

Since the formulae (l55|) - ( !56l) are very complicated, we organize the calculations by the 
diagrammatic method. Since Feynman diagrams have not been used before in inverse spec- 
tral theory, we digress to present the fundamentals of the diagrammatic approach to the 
stationary phase expansion; clear expositions are given in [Al [E] (see also [AG] ) . 



5.1. Stationary phase diagrammatics. We consider a general oscillatory integral 



where a G C^(]R") and where S has a unique critical point in supp a at 0. We write H for 
the Hessian of 5* at and for the third order remainder in its Taylor expansion at x = 0: 

S{x) = S{0) + {Hx, x)/2 + Rsix). 

The stationary phase expansion is: 

- 1^1=0 z^v=o["'ydj)\-v\y^3[-gj)) \j=o jj • 
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The graphical analysis of the stationary phase expansion consists in the identity 

(^^^ ^''^W^W^'^W^ — 7^ — - ^ \Mm 

where Gyj is the class of labelled graphs {Q-,t) with V closed vertices of valency > 3 (each 
corresponding to the phase), with one open vertex (corresponding to the amplitude), and 
with / edges. The function i 'labels' each end of each edge of Q with an index j G {1, . . . , n}. 

Remark 5.2. The term ' open vertex' is equivalent to ' marked' or 'external' vertex in 
some texts, and is graphed here as an unshaded circle. A ' closed' vertex is the same as an 
'unmarked' or 'internal' vertex and is graphed as a shaded circle. Also, it is non-standard to 
include the labels i in the notation for Feynman amplitudes; we do so because in our problems 
certain labels are distinguished. 

Above, |ylMt(^)| denotes the order of the automorphism group of Q, and Ie{Q) denotes the 
'Feynman amplitude' associated to the labelled graph By definition, Iii^Q) is obtained 

by the following rule: To each edge with end labels m, n one assigns a factor of -^h"^"" where 

as above H^^ = (/i*""). To each closed vertex one assigns a factor of ^^^^r^^^ir where v is 
the valency of the vertex and ii . . . ,1,^ at the index labels of the edge ends incident on the 
vertex. To the open vertex, one assigns the factor g^fi " al'.^ ' where u is its valence. Then 
is the product of all these factors. To the empty graph one assigns the amplitude 1. 
In summing over {Q, €) with a fixed graph ^, one sums the product of all the factors as the 
indices run over {1, . . . , n}. 

We note that the power of in a given term with V vertices and / edges equals k'^s' ^ 
where Xg' = V — I equals the Euler characteristic of the graph Q' defined to be Q minus the 
open vertex. We thus have; 

We note that there are only finitely many graphs for each x because the valency condition 
forces / > 3/2V. Thus, V < 2j, I < 3j. 

5.1.1. Stationary phase formula for I^}^^ . Since Feynman diagrams and amplitudes are un- 
familiar in wave trace calculations, we digress to give some details of the proof of (1^51) and 
to tie it together with the form of the stationary phase expansion in standard texts in par- 
tial differential equations (cf. [Ho] I). This latter form can also be used to corroborate the 
calculations below. 

The stationary phase of ( |Hoj I. Theorem 7.7.5) reads: 

(55) z^-i'-^r'^ ' j:k-'VAo) 



y^\detH\ 



j=0 



where 



(56) ^.a(0)=E E n77^(^"'^'^)'^W)U=o 

u—p=j 2u>3fj. ^' 
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In diagrammatic terms, the pair (/i, u) correspond to graphs with u = I edges and fi = 
V closed vertices, hence of Euler characteristic fi — u = —j. We note that the factor 
is common to all graphs of Euler characteristic —j and in our analysis we absorb into 
the prefactor. To tie (156|) together with (153|1 . we sketch the proof of the latter, following 
the exposition in [E] in the case where the amplitude is = 1. We outline the procedure 
following the notes of Etingof [Ej This special case turns out to be the most important 
for the applications in this paper, since terms with derivatives of the amplitude will not 
contribute to the highest order jets in the wave invariants. The notes of Axelrod [AJ give a 
clear discussion (as above) of the contribution of the amplitude to the Feynman amplitude. 

Proposition 5.3. We have: 



■{H-'D,Dr{Ri;)U=o= E 



fi\u\' ' ' ^^■'^^1-='^ ^ \Aut{g)y 

Proof. We need to re-write the left side as a sum over graphs in Gu,^i (the class of graphs 
with edges, /x closed vertices of valency > 3). 

Let n = (no,ni,...) be a sequence of non-negative integers, of which all but a finite 
number are zero, and let G'(n) denote the set of graphs with tiq 0-valent vertices, ni 1-valent 
vertices etc. We are only considering the case where the amplitude equals one, one so there 
are no external vertices. 

We write Rsi^x) = X]m>3-^m(^' • • • )3;)/m!, where Bm = d"^S{0), as a sum of its homo- 
geneous terms. Change variables x —>■ Vkx, write e'^^^^Vk^^ = Ji^e^^""^^^^^^' and Taylor 
expand each exponential to obtain 

= En ^n, with 

(57) 

The integral may be calculated by Wick's formula. The diagrammatric interpretation at- 
taches to each factor iBm a 'flower' of valency m, i.e. a closed vertex with m outgoing edges. 
Thus, the index n prescribes a set of Um flowers of valency m. Let T be the set of the ends 
of the outgoing edges of all of the flowers. For each pairing a of the ends one obtains a graph 

Associated to each graph is its Feynman amplitude -Fn o-- As described above, one labels 
each end of each edge of the graph by indices in {1, . . . , n}, assigns a factor of -^hJ^^ to an 
edge with end labels m, n and flower (closed vertex) of valency i with end labels (a;„^, . . . , a;„.) 
one assigns a factor of ife g^f^'^^g^n. . One multiplies these expressions over all edges and closed 
vertices and then sums over all labelings. One then has 

_ (27r)"/^ I (rr.^^^ ^ 

VdetH (m!)"-n„! ^ 
By comparison, in (l56l) . one Taylor expands the full factor e^^ to obtain 
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Since 

(58) 1 V (^]u V n ^ 

n:ln|=/i n:|n|=/j 

it follows that 



m 



(59) ^{H-'D,Dr{R^,)U=o = ^-^{H-'D,Dr Yl ^ ^ 



LL\v\ ' . - ^, . . ,i — ^ (m!)"'"(nm)! 

n:|n|=/^ 

For each fixed n, the term on the right side for this n is the i^th term in the expansion of 
Zn when (as in the proof in [Ho] ) one applies the Plancherel formula to the integral (157|) for 
and Taylor expands e*^ hj^y)/"^ ^ The z/th term can be sifted out by replacing H \H and 
finding the term of order X~'^ on each side. Note that [fi, v) are determined by n : Indeed, 
/i = Ylim ^rn, and since each outgoing vertex is paired with exactly one other outgoing vertex 
to form an edge, = | ^m^^m- We write /i(n), i^(n) for the these values. The terms 
in the sum over n with |n| = fi run over those n for which i^(n) = z/, and thus we have 

-{H ^D,DY J2 ^^frnl)-^(nrny. ^^""(mir-nj. ^ 

n;|n|=/x n:|n|=/j,!^(n)=!^,(7 

Finally, as explained in 

The same identity holds if we restrict to pairings and graphs with /i vertices and u edges. 
Cancelling common factors, we get 

Combining with fl59l) completes the proof. 

□ 

5.2. Maximal derivative terms. We now apply the diagrammatic stationary phase method 
to the oscillatory integrals f|T9|) . Further, we consider the additional aspect of extract- 

ing from the stationary phase expansion the terms which involve the highest number of 
derivatives of the boundary defining function /± in each power of k~^. Such terms with 
the maximal number of derivatives arise only from special graphs and from special terms 
in the corresponding Feynman amplitudes with special labelings of the vertices. This is a 
non-standard feature of diagrammatic analysis and indeed depends on the very special phase 
and amplitudes in I"^'^^ ■ A further key issue is the dependence on the number of iterates M 
of the bouncing ball orbit. 

For emphasis, we state our objective as follows: 

• Enumerate the diagrams of each Euler characteristic whose amplitudes contain the 
maximum number of derivatives of dfl among diagrams of the same Euler characteris- 
tic. Determine which vertex labellings produce the maximum number of derivatives. 
Then determine the corresponding "maximal derivative Feynman amplitudes", i.e. 
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the sums of monomials containing the highest number of derivatives. We denote 
them by /™^^(^). 

As we will see, only the principal oscillatory integrals of Definition 14.31 give rise to terms 
in I'^^^IQ). We use the following notation for the class of labelled graphs which give rise to 
two types of maximal derivative terms. 

• ^ul'fl'^ *^ the (not necesssarily unique) labelled graphs whose Feynman ampli- 
tude contains terms of the form /'•'^•*(0)/*^*-'(0)aQ'^^(0). In fact, we will show that c = 
for all labelled graphs contributing to the highest number of derivatives of / in a 
given order of wave invariant. 

We denote by J"^ the operation of extracting the terms with p derivatives. That is, J'^ 
applied to a monomial in derivatives of the phase is equal to the monomial if it contains 
a factor with p derivatives of the phase and zero otherwise. From Proposition 15.31 we can 
evaluate the combinatorial coefficients of Feynman amplitudes with a specified number of 
derivatives. 

Corollary 5.4. We have: 

^''a^^^^^^)Uo=xi=-=x,™=o= \Aut(g)\- 

5.3. The principal terms. Our first step is to analyze the stationary phase expansions of 
the principal terms l2r'p'^i^) sense of Definition 14.31 By Proposition 14.41 it suffices to 

consider We show that the non-principal terms only contribute lower order derivative 
data to the Balian-Bloch invariants B^j. In the next section, this data will be proved redun- 
dant in the case of the symmetric domains of this article. As mentioned in the introduction, 
we only use the attributes of the phase and amplitude described in Theorem 14. 2[ We now use 
this information to determine where the data f±{0), f^''~^\o) first appears in the stationary 
phase expansion for the oscillatory integrals. 

The only critical point occurs where x = 0. We denote by 7i± the Hessian operator 
in the variables {xi, . . . ,X2r) at the critical point x = of the phase C±. That is H± = 
{Hess{/:±)-^D, D), where D is short for ■■■ , ^J-). 

5.3.1. The principal term: The data f±{0). We first claim that f^''\o) appears first in the 
k~^^^ term in the stationary phase expansion of ^r'p'^ ■ This is because any labelled graph 

{Q, i) for which Ie{Q) contains the factor f^''\o) must have a closed vertex of valency > 2j, or 
the open vertex must have valency > 2j — l. The minimal absolute Euler characteristic 
in the first case is j — 1. Since the Euler characteristic is calculated after the open vertex is 
removed, the minimal absolute Euler characteristic in the second case is j (there must be at 
least j edges.) Hence such graphs do not have minimal absolute Euler characteristic. More 
precisely, we have: 

Lemma 5.5. In the stationary phase expansion of 12°'^^^ , the only labelled graph {Q,i) with 
—x{0') = J — 1 with Ie{G) containing f^''\o) is given by: 

• Gi'^f'^ G Gij (i.e. ^ = V = l,I = v = i). There is a unique graph in this class. It 
has no open vertex, one closed vertex and j loops at the closed vertex. 
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Figure 5. x = j — 1,V = 1,1 = j), j loops at one closed vertex. All 

labels the same. Form of Feynman amplitude: {W^y D^xp Cj^ = {h^y /^"^^^O) 

• The only labels producing the desired data are those ip which assign all endpoints of 
all edges labelled the same index p. 
The J'^Hh part of the Feynmann amplitude is 

I-^-i^f) = ArL iw^,) A(r){ {h'^Y ff\Q) - {hl'V f^''\o)} , 

where we neglect terms with < 2j — 1 derivatives. 

We are also interested in the f^'^~^\o) terms, but postpone the calculation of the f^''~^\o) 
- terms arising from the diagram until Lemma [5.6( ii) (they turn out to vanish). 

Proof. By ( !56|) . the data f±{0) only occurs in the term = 1, z/ = j of (!56|) . To see this, we 
note that the Hessian operator TY*^ associated to has the form 



dxi-^ Oxj-^ . . . Oxi^ dx 



Any term {h^Dl^Y applied to produces a f^^\o) term. 

We can also argue non-diagrammatically that no Uj > 2{j + 1), i.e. the power k~^^^ is the 

greatest power of k in which f^^\Q) appears. Indeed, it requires 3/i derivatives to remove 
the zero of i?3. That leaves 2z/ — 3;U = 2j — 2 — /i further derivatives to act on one of the 
terms D^R^, or 2j — 2 — derivatives to act on the amplitude. The only possible solutions 
of (z/, /i) are (j — 1,0), (j, 1). Referring to statement (i) of Theorem 14.21 and to we see 
that the principal symbol of the amplitude depends only on /±,/4, so there is no way to 

differentiate the amplitude 2j — 2 times to produce the datum f^''\o). Hence, (z/, fi) = (j, 1) 
and the only possibility of producing f^-'\o) is to throw all 2j derivatives on the phase. 

Now let us determine I^^^^^iQ) for the labelled graphs {Q, i) above. The terms with maximal 
number 2j of derivatives in the Feynman amplitude (apart from the overall universal factor 
in dH])) are given for some non-zero constant Cg by 

(61) 

= Cg{4rL)Ao{r) ElUm'^^MftliO). 

The factor {4rL)Ao{r) comes from the leading value of the amplitude (cf. Lemma [4. 5p . By 
Proposition [531 Cg = jj^^^ = wg. 
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Indeed, to obtain /± (0), all labels at all endpoints of all edges must be the same index, 
or otherwise put only the 'diagonal terms' of Ti'^, i.e. those involving only derivatives in a 
single variable can produce the factor f^-'\o). We then use Lemma (va) to complete 

the evaluation. The part of the pth term {h^Y Dl^^C+{0) of the sum which involves /^^^(p)(0) 
equals 

by (|49l). 

We then break up the sums over p of even/odd parity and use Proposition 12.51 to replace 
the odd parity Hessian elements by and the even ones by h'^. Taking into account that 
w+{p) = 1(— 1) if p is odd (even), we conclude that 

(62) ^ 8rL (u;^.„o) Ao{r){ {h^y f^^'\o) - {h'y f^^'\o)} + ■■■, 

where again ■ ■ ■ refers to terms with < 2j — 1 derivatives. We observe that, as claimed, 
the result is invariant under the up-down symmetry /_|_ <^==^ — /_ and under the right left 
symmetry /±(x) ^ f±i-x). 

a 

Thus, we have obtained the even derivative terms in Theorem 15.11 

5.3.2. The principal term: The data ^\^)- We now consider the trickier odd-derivative 
data f'^^ ^''(0) in the stationary expansion of l2r'p^, which will require the attributes of the 
amplitude (l45l) detailed in Theorem 14.21 

We again claim that the Taylor coefficients f^'' ^\o) appear first in the term of order 
/^-i+i Further, only five graphs can produce such a factor, and of these only two contribute 
a non-zero Feynman amplitude. These two graphs are illustrated in the figures. In the 
following section, we will show that /jf"* ^''(0) can only occur in higher order terms in 
also in the singular trace terms. 

To prove this, we first enumerate the labelled graphs Q in the stationary phase expansion 
of /2rp"* whose Feynman amplitude Ie{G) contains a factor of f^''~^\o) in the term of order 
k~^~^^, and we show that this data does not appear in terms of lower order in k~^. 

We recall that = means equality modulo -R2r(i7^-'~^/+(0), jr^-'~^/-(0)). 

Lemma 5.6. In the stationary phase expansion of l2°'^^ , 

(i) There are no labelled graphs Q with — x'(^) := —x{G') < J ~ 1 for which h{Q) contains 
the factor f^^^-^\0). 

(ii) There are exactly two types of labelled diagrams with x{G') = — J + 1 such that 
h{Q) is non-zero and contains the factor f^''~^\o). They are given by (see figures): 

• ^2^+1^'" ^ ^2j+i with V = 2, 1 = j + 1: Two closed vertices, j — 1 loops at one closed 
vertex, 1 loop at the second closed vertex, one edge between the closed vertices; no 
open vertex. Labels ip^g: All labels at the closed vertex with valency 2j — 1 must be the 
same index p and all at the second closed vertex must the be same index q. Form of 
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Thus, this graph contributes 

p,q=l 

• Q2'!j+i'^'^ ^ with V = 2, 1 = j + 1: Two closed vertices, with j — 2 loops 
at one closed vertex, and with three edges between the two closed vertices; no open 
vertex. Labels ip^g: All labels at the closed vertex with valency 2j — 1 must be 
the same index p and all at the second closed vertex must the be same index q; 
{hfy-\hP^fDli~'C±Dl^C± = {hP^y-\h^^Yf^^'-'\0)fi\0). Thus, this graph con- 
tributes 

2r 
p,q=l 

• In addition, there are three other graphs whose Feynman amplitudes contain factors 
off^^~^\0). But for our special phase and amplitude, the corresponding amplitudes 
vanish. 

Proof. It will be seen in the course of the proof that only connected graphs can contribute 
highest order derivative data (the amplitude for a disconnected graph is the product of the 
amplitudes over its components). Connected labelled graphs with — x' < j — 1 for 

which Ie{Q) contains the factor f^'' ^\o) as a factor must satisfy the following constraints: 

• (a) Q must contain a distinguished vertex (either open or closed). If it is closed it 
must have valency > 2j — 1. If it is open, it must have valency 2j — 2. We denote by 
i the number of loops at this vertex and by e the number of non-loop edges at this 
vertex. 

• (b) -x{Q')=i-v<j-i 

• (c) Every closed vertex has valency > 3; hence 21 > 3V. 

We distinguish two overall classes of graphs: those for which the distinguished vertex 
is open and those for which it is closed. Statement (a) follows from the attributes of the 
amplitude in Theorem 14.21 In the first case, 2j — 2 derivatives must fall on the amplitude 
(i.e. the open vertex) to produce f^^~^\o). In the second case, 2j — 1 derivatives must fall 
on the phase (i.e. the closed vertex). 

We first claim that V < 2 under constraints (a) - (c). When the distinguished vertex is 
open, then \^ = if — x' = J — 1 (as noted above), and there are no possible graphs with 
—x' < J ~ 2. So assume the distinguished vertex is closed. Let us consider the 'distinguished 
flower' To consisting just of this vertex and of the edges incident on it. Denoting the number 
of loops in To by i, we must have 2£ + e > 2j — 1 edges in Tq to produce f^-'~^\o). We then 
complete Fq to a connected graph Q with — x' < j — 1. We may add one open vertex, V — 1 
closed vertices and N new edges. 
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Suppose that there is no open vertex. We then have: 

(i) 2£ + e > 2j - 1 

(63) i {ii) i + e-V + N = ] - I 

(m) e + 2N >?,{V - 1) 

The last inequahty follows from the facts that each new vertex has valency at least three, 
and that each of the r edges begins at the distinguished vertex. Solving for V in (ii) and 
plugging into (iii) we obtain < 3j — 3£ — 2e. Plugging back into (ii) we obtain V < 
2j - 2£ - e + 1 < 2j + 1 - (2j - 1) = 2, by (i). Thus the claim is proved. 

Now suppose that Q contains one open vertex and V closed vertices. Then (i) and (ii) 
remain the same since the xiG') is computed without counting the open vertex. On the 
other hand, (iii) becomes e + 2N > 3{V — 1) + 1, since the open vertex has valence at least 
one. This simply subtracts one from the previous computation, giving V < 1. Thus, the 
distinguished vertex is the only closed vertex. 

Now we bound in the connected component of the distinguished constellation. First 
suppose that V = 1. There is nothing to bound unless the graph also contains one open 
vertex, in which case N counts the number of loops at the open vertex. We claim that = 
in this case. Indeed, we have i + e + N = j . Substituting in (i), we obtain 2A^ + e < 1. The 
only solution is = 0, e = 1. 

Next we consider the case V = 2. As we have just seen, no open vertex occurs. From (i) 
+ (ii) we obtain 2A^ + e < 3, hence the only solutions are A^ = 1, e = 1 or A^ = 0, e = 3. 

We tabulate these results as follows: 



Graph parameters 


V 


i 


e 


N 
















1 


1 


j 











1 


j-1 


1 





1 


2 


j-1 


1 


1 





2 


j-2 


3 









We now determine the Feynman amplitudes for each of the associated graphs. As we will 
see, the amplitudes vanish for the first three lines of the table, and do not vanish for the 
last two. The non- vanishing diagrams are pictured in the figures (Figures 6 and 7). 

• (i) The only possible graph with V = is: Qq^^J^^, V = 0, 1 = j — 1: j — 1 loops at 
the open vertex. Taking into account the structure of the amplitude in Theorem 14. 2[ 
in order to produce /^^•'"^^(O), all labels at the open vertex must be the same index 
p. We claim that the Feynman amplitude vanishes: 

2r 

(64) = (Const.) Y,{hn'-'Dt'A ^ x ft'''\o) = 0. 
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Indeed, this is the case (/x, z/) = (j — 1, 0) of ( l56l) . which corresponds to applying all 
derivatives Dl^~^ on the principal symbol aP of the amplitude for some p = 1, . . . , 2r, 
and it is proved in §4.11 flSOl) that it vanishes. 

• (ii) QI^J' C Qi,j,V = 1,1 = j: j loops at the closed vertex. This is the graph 
which produced /'■^-'^(O), and we now verify that it does not produce an amplitude 
containing f^^-'~'^\0). To produce /'•^-'"^-'(O), all but one label must be the same (p), 
the last label different {q ^ p)- Feynman amplitude: 

2r 

™ = {Const.) V {wy-^wm^^^-^^D^^c = {w^y-^w"'f^2'-^\<S)f^{Q) = o. 

p,q=l 

The vanishing is verified in §4.11 flSTj) . 

• (iii) C Qij, V = 1, 1 = j: j — 1 loops at the closed vertex, one edge between 
the open and closed vertex. To produce /'•^-'"^''(O), all labels at the closed vertex 
must be the same index p. We claim that again the Feynman amplitude vanishes: 

2r 

I^,fL,_, = {Const.) {hPP)^-^hPW^J~^CDy = x /('^"'^(O) = 0. 
p,(?=i 

Indeed, exactly one derivative is thrown on the amplitude. To check this, we note 
that this is the case (/i, i^) = (j, 1) of ( l56i) in which is applied to a^R^. To 
produce the data f^'' ^\o), the operators Dl?~^Dx^ contribute by applying D^^^^ to 
i?3 {p = 1, . . . , 2r), and by applying the final derivative D^^ to the amplitude. But 
VflO (0) = by dSD- 

• (i'^) ^I^Z+f'^ <^ ^2,j+i(— X = j — 1; ^ = 2, / = j + 1): Two closed vertices, j — 1 loops 
at one closed vertex, 1 loop at the second closed vertex, one edge between the closed 
vertices; the open vertex has valency 0. All labels at the closed vertex with valency 
2j — 1 must be the same index p and all at the closed vertex must the be same index 
q. Since there are no derivatives of the amplitude, we extract its principal term and 
obtain 

= 8rLAr{0)C^..-^Y''^7^^^^^^^^ 

The calculation of the coefficients is similar to that in (iii), except that now we have 
two factors of the phase. The factor containing 2j — 1 derivatives of C is evaluated 
in (iv) - (v) of the table in Lemma 14.51 and the third derivative factor is evaluated in 
§4.1.11 (4). Again the combinatorial constant is evaluated in Proposition 15. 3[ 

• (v) There is a second graph C ^2j+i(— X = j ~ 1; ^ = 2, / = j + 1): It has 
two closed vertices, with j — 2 loops at one closed vertex, and three edges between 
the two closed vertices; the open vertex has valency 0. Labels ip^q-. All labels at the 
closed vertex with valency 2j — 1 must be the same index p and all at the closed 
vertex must the be same index q. Again, there are no derivatives on the amplitude. 
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and we get 



max 



g. 



(65) 



^ 2rLAMC^-.., E?,=i(/^T)^'-^(/^^+^) V(j>)«^+(g)£Sj^(0) 

As noted above (cf. §4.11 (!52l) ). other (mixed) third derivatives of C vanish on the 
critical set. The combinatorial constant is evaluated in Proposition 15.31 

We now combine the terms in (iv) and (v) and evaluate the coefficients to obtain 

2rLA(0) {(«;,..-M.) Y.l,=,my-'h':h'-^ 

We obtain the expression stated in Theorem (15.11) by breaking up into indices of like 
parity and using Proposition 12.51 

□ 



We pause to review the sources of the various constants and to check that sums over the 
several ± signs do not cancel. In particular, it is crucial that the coefficient of 

non-zero, since it is this term which determines odd Taylor coefficients and allows us to 
decouple even and odd derivative terms. 

Remark 5.7. The constants and sums over ± are of the following kinds: 

• The factor of C in the amplitude produces 2rL. 

• The following ± signs arise (with some redundancy): 7^, f±, w± or equivalently 
C±, p even (odd), and the two terms of C which depend on a given index Xp (J^- 
Proposition \4.4\ shows that the two possible choices of w± produce the same data. 
Since 7 = 7""'^ there is no question of cancellation between B^± . 

• The odd derivative monomials with maximal derivatives of f have the form 

^(o)/f (0), ff~'\o)ff\o), /i'^-^^(o)/f (0), f'^'-'\o)ff\o). 

By Theorem \ 5.1\ the wave invariants are invariant under —f-,f- — ^ ~/+; 

hence the only possible cancellation could occur between /|''"'^(0)/f (0) and /|'^"'^(0)/i' 
However, no such cancellation occurs, as noted after the calculation in [4^ )> or in 
Theorem \5.1\ where it is noted that the monomials always occur in the form 

In fact, the ± sum in each factor D'^-'^C, Dl?~^C, D^^C gives rise to a factors of 4 in 
odd derivative terms, and factors of 2 in even derivative terms. For the same reason, 
no cancellations occur between the sum over p even versus p odd. 
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Figure 6. 



^2j-l,3,0 



2,i+l 



-X = j-l;V = 2,I = j + l] 



Figure 7. {v) : g^^r^Y'' c e;2,,+i(-x = j - 1; 1^ = 2, / = j + 1) 

5.4. Non-principal terms. To complete the proof of Theorems 14.21 and 15.11 it suffices to 
show the non-principal oscillatory integrals J^}'^ with M > 2r do not contribute the data 

f^''\o), ^^(0) to the coefficient of the k~^^^- term (or to the k~"^ term for any m < j — 1). 

We recall from Proposition 13.101 that /^'^ can only have a critical point if M > 2r and 
M — \a\ = 2r. In the non-principal terms where M > 2r, the oscillatory integral /^'^ is 
obtained by regularizing the kernel of in Proposition 13.81 which is an oscillatory integral 
with a singular phase and amplitude (cf. [Z5], §6). 

The regularization produces the oscillatory described in Corollary 13.91 In the case where 
M — |(t| = 2r it is an integral over T^^ with the same phase as in the principal terms but with 
an amplitude of order — |cr|. The sum over M in Proposition 13.61 and over a in (1291) can thus 
be seen as the construction of an oscillatory integral expression for the trace of Proposition 
13. 6[ with an amplitude obtained by regularizing the sum of singular oscillatory integrals. 

The stationary phase analysis of the sub-principal terms J^}'^ is therefore almost essentially 
the same as for the principal term. The only additional feature is the following description 
of the amplitude: 
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Lemma 5.8. The amplitude Afj{k, ipi, ip2) of N^r in Proposition 3.8 is a semi- classical ampli- 
tude of order —\a\. In its semi-classical expansion A„{k,Lpi,Lp2) ~ Xl^o ^ ''^'""^^."■(v'l' V's); 
the term A^^n depends at most n + 2 derivatives of f . In particular, the value D^Aa^n\^c 
of its ath derivative at the critical point depends at most on n + 2 + |a| derivatives of f at 
X = 0. 

Proof. The algorithm for calculating y4o-(fc, V52) is given in [Z5] §6 (see also [AGj ) . We 
briefly review the algorithm in order to prove that the amplitude has the stated properties. 

The algorithm consists in successively removing factors of Nq from compositions of Nq 
and A^^i in N„ (cf. §3). The first step consists in expressing the compositions Nq o A^^i and 
A^^i o A^Q as oscillatory integrals of one lower order (cf. Lemma 6.2 of |Z5j ). From the explicit 
formula for the composition (cf. (74) of [Z5j). the new amplitude A{k + iT,(pi,(p2) has the 
form 

(66) A{k + ir, (pi, (p2) = / xik,u, (pi, (p2)G{k + iT,u, (fi, (p2)\u\Hi^\{k + iT)\u\)e''"'''du, 

Jr 

where x is a suitable cutoff and G is a semi-classical amplitude constructed from the ampli- 
tude of Ni (cf. (78)-(79) of [Z5]). Also, a = sin((g((/?2) - qM, i^gi^^))- 

The amplitude G is constructed as follows: From A"o one obtains a contribution of 
H[^\{kfi + iT)\q{ip3) — q{ipi)\) cos Z{q{ip3) — q{(pi), 1^9(^3)), while from A"! one obtains a semi- 
classical amplitude. One changes variables by putting 



(67) 




|g((^3) - g((^i)|, > <^3 



under which the amplitude of A^i is transformed to a smooth amplitude of the same order 
in {ip2,u), while the factor of cos Z{q{ips) — q{fi), ^q('P3)) changes to \u\K{(pi,u) where K is 
smooth in u. A simple calculation shows that K{ipi,0) = —^n{ipi). The full amplitude G 
is a product of these two factors. One sees that it depends analytically on /, /', /" with /" 
coming from the cosine factor. 

One then Taylor expands G in m and verifies that it produces a semi-classical expansion of 
A{k + ir, ipi, ip2)- The du integrals can be explicitly evaluated using the cosine transform of 
the Hankel function ( |Z5] . Proposition 4.7; see also [AGj ) . The \u\du in the cosine transform 
gives rise to a factor of and the factor of A'o carries a factor of k, so that the removal of 
Nq introduces a net factor of k~^. This factor is responsible for the lowering of the order by 
one for each removal of Aq. 

The coefficient of A;~^~" in the final amplitude thus derives from the nth term in the Taylor 
expansion of G{k,u,ip) in u and in particular depends on the same number of derivatives 
of /. Since G is an analytic function of /, /', /", it follows that the k~^~^ term depends at 
most on n + 2 derivatives of /. 

The process then repeats as another factor of A"o is removed from the resulting composition. 
The same argument shows that each elimination of A"o introduces a new factor of k~^ which is 
unrelated to Taylor expansions of G. We now verify that after r repetitions of the algorithm, 
the new amplitude is semi-classical and its k~^'~^ term depends on only n + 2 derivatives of 
/■ 
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We argue by induction, the case r = 1 having been checked above. After r — 1 steps, 
we obtain an oscillatory integral operator with an amplitude Ar-i satisfying the hypothesis 
and with the phase of A^i. We then apply the algorithm for the composition of Nq with this 
oscillatory integral operator. It has the form of (166!) except that now G = Gr is constructed 
using Ar-i and Nq. The algorithm is to multiply A^—i by the cosine factor above, to change 
variables to u, to Taylor expand the cosine factor to one order to obtain \u\K and to define 
Gr = KAr-iJ where J is the Jacobian. The Taylor expansion producing K is responsible 
for the initial increase in the number of derivatives of / to /". After that point, it is only 
the Taylor expansion of G^ in u which produces further derivatives of /. Thus, the number 
of derivatives of / in the term of order is n + 2. 

It follows that, after removing all \a\ factors of Aq, one obtains an amplitude which is of 
order — |cr| and whose term involves at most n derivatives of /". 

□ 

Lemma 5.9. The non-principal terms do not contribute the data f±{0), f'^~'^{0) to the term 
of order k"^"^ . 

Proof. We consider the diagrammatic analysis of J^}'^ along the same lines as for the principal 
term. The only new aspect is the amplitude. Since it now has order —\a\ < 0, the terms 
where one differentiates the phase to the maximal degree now have order k^i+^^\°'\ and thus 
do not occur in the k~^~^ term. 

The only remaining possibility is that the data could occur in terms where one differentiates 
the amplitude to the maximal degree. By Proposition 15. 8[ the term of order fc"!*^!"" contains 
at most n + 2 derivatives of /. To obtain a term of order —j + 1, one needs |cr| + n < j — 1 
and one can take only 2(j — 1 — |cr| — n) further derivatives in the k^^^^ term. This produces 
a maximum of 2j — 2|cr — n derivatives of /. The maximum occurs when n = 0, in which 
case there are < 2j — 2|o"| < 2j — 2 derivatives of /. 

□ 

For emphasis, we determine the lowest order term in which such data do occur: 

SuBLEMMA 5.10. In the stationary phase expansion of the non-principal term I^f^, the data 
f± {0), f±~^{0) appear first in the k^~^~^'^^- term. 

Proof. To determine the power of k^^ in which this data first appears, we need to minimize 
|cr| + z/ — /i subject to the constraint that 2i/ — 3/i > 2j — 3. This is \a\ plus the constrained 
minimum of u — fi. The sole change to the principal case is that the constraint is 2z/ — S/z > 
2j — 3 in the top order term of the amplitude rather than 2z/ — 3/i > 2j — 2. Since the 
solutions must be non-negative integers, it is easy to check that again ly > j — 1 and that 
(yLi, z/) = (0,j — l),(l,j) achieve the minimum of u — fi = j — 1. If there are r drops 
in the symbol order, we need to minimize \a\ + r + i/ — /i subject to the constraint that 
2z/ — 3/i > 2j — 3 — r. The minimizer produces the result stated in the Sublemma. 

□ 

This completes the proof of Theorems 14.21 and (15.11) . 

5.5. Appendix: Non-contributing diagrams. In figures (6)-(7), we displayed the dia- 
grams which contribute non-zero amplitudes to the leading order derivative terms. For the 
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Figure 9. (ii) : c Gij). 

sake of completeness, we also include diagrams which do not contribute because the corre- 
sponding amplitudes vanish. Figures (8) - (9) are labelled consistently with the discussion 
above. Figure (5) is also a 'non-contributing diagram' to the leading order odd derivative 
term. 
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5.6. Balian-Bloch invariants at bouncing ball orbits of up-down symmetric do- 
mains. We now simplify the expression in Theorem 15. II in the case of Z2-symmetric domains. 
The following result, stated in ([9]), is essentially a corollary of Theorem 15. 1[ It uses one sim- 
plification which will be proved in Proposition 16. 5[ 

Corollary 5.11. Suppose that (fi, 7) is invariant under an isometric involution a, and 
that 'J is a periodic 2-link reflecting ray which is reversed by a. Then, modulo the error term 
-R2r(i7^"'~^/(0)), Byj-i is given by the expression 

Proof. Using that /_ = —/+, we can cancel the signs in the formula of Theorem 15.11 and add 
the top and bottom to obtain, 

^ ^tL Ao{r) {{Wg.^,o) EUih^' f^^'KO) 

+ 4(^^;^2.-l,3,o)(/^Pf)^-2(/^f9)3]}/(3)(0)/(2J-l)(0). 

Further, in this Z2-symmetric case, all of the coefficients h^^ are clearly equal. The sum 
Yl^Li h^'^ is independent of p and is evaluated in Proposition (16. 5p . leaving the stated ex- 
pression. □ 

6. Proof of Theorem ( 11.11) 

We now prove the inverse spectral result for simply connected analytic plane domains with 
one special symmetry that reverses the endpoints of a bouncing ball orbit. The method is 
to recover the Taylor coefficients of the boundary defining function from the Balian-Bloch 
invariants at this orbit. 

As simple warm-up for the proof, we give a new proof that centrally symmetric convex 
analytic domains whose shortest orbit is the unique orbit of its length (up to time-reversal) 
are spectrally determined within that class: 

Proof of Corollary 11.21 Consider the wave invariants of the shortest orbit as given in 
Theorem 15.11 They are spectral invariants since the shortest length is a spectral invariant. 
By Ghomi's theorem [Ghj . the shortest orbit is a bouncing ball orbit. The orbit must be 
invariant under the two symmetries up to time-reversal since its length is of multiplicity one. 
Hence, the two symmetries imply that /+ = — /_ := / and that /'^^■'"'"^^(O) = for all j. It 
follows that /^^"'^(O) are spectral invariants for each j, and thus the domain is determined. 
QED 

The same proof shows that simply connected analytic domains with the symmetry of an 
ellipse and with one axis of prescribed length L are spectrally determined in that class. 

6.1. Completion of the proof of Theorem II. li We now complete the proof of Theorem 
11.11 Thus, we assume that (fi, 7) is up-down symmetric, i.e. is invariant under an isometric 
involution a, and that 7 is a periodic 2-link refiecting ray which is reversed by a. 

There are two overall steps in the proof. First, and foremost, we study the expressions 
in Corollary 15.111 The key point is that the Hessian of the length function is a circulant 
matrix in the symmetric case, and that allows us to analyze the Hessian sums which occur 
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as coefficients in the Balian-Bloch wave invariants. In particular, we decouple even and odd 
derivatives using the behavior of the Hessian sums under iterates 7*". After that, a simple 
inductive argument shows that all Taylor coefficients of /+ may be determined from the 
Balian-Bloch invariants. 

We now begin the analysis of the Hessian sums. 



6.2. Circulant Hessian at Z2-symmetric bouncing ball orbits. In the case of Z2- 

symmetric domains in the sense of Theorem ll.il Ra = Rb '■= R and 

(68) cosa;/2 = 2(l — — ) (elliptic case), cosha/2 = 2(l — — ) (hyperbolic case). 

R R 



We put: 
(69) 



-2cosa;/2 (elliptic case), a = — 2cosha;/2 (hyperbohc case) 



By [16] and Proposition [2?2l the Hessian of the Length function in Cartesian graph coordinates 
simplifies to: 



(70) 



a 1 



1 a 







H2r = < 



1 a 1 
1a 1... 

^1 ... a 

We observe that fITOl) is a symmetric circulant matrix (or simply circulant) of the form 
(71) (-L)/72r = C(a,l,0,...,0,l), 

where a circulant is a matrix of the form (cf. [D] ) 



( Cl C2 



(72) 



Cfi Cl ... Cn~\ 



C(ci,C2, . . . ,C„) = < 



> . 



C2 C3 



Cl 
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Circulants are diagonalized by the finite Fourier matrix F of rank n defined by 

r 1 1 ... 1 



(73) 



p* 



n 



-1/2 



W 



w 



n-l 



I 1 



W 



n-l 



W 



(„_!)(„_!) 



W 



Here, F* = (F)^ = F is the adjoint of F. By [D], Theorem 3.2.2, we have C = F*AF where 



(74) A = Ac = diag (pc(l), • • . ,Pc(«^"~ ), with pdz) = ci + + ■ ■ ■ + c^z 
Here, by diag we mean the diagonal matrix with the exhibited entries. 

6.3. Diagonalizing H2j}. Applying the above to C = H2r- 
Proposition 6.1. We have: 



n-l 



H2r 



-L F*{dtag{^^,... 



1 , r, (2r-l) 

a+2 cos '- 



l)7rs )-^) 



where a is defined in f6W . 

Proof. We use the notation Pa,r{z) for Pc{z) in the case where C = C(a, 1, 0, 



0, 1). Thus, 



(75) 

By fl73|) we have, 
(76) 

Since 
(77) 

we have 
(78) 

and inverting 



Pc[z) 



a + z + z 



2r-l 



H. 



2r 



FMiag (Pa,r(l), • 



,Pa,r(W 



2r-l 



[w 



e >■ 



Pa Aw 



a + w+w~ 



[w 



e ^ 



H2r 

nves the statement. 



— 1 {2r — IW 

— FMiag (a + 2, . . . , a + 2 cos ^ '—)F, 

L r 



□ 



6.4. Matrix elements of at a Z2-symmetric bouncing ball orbit. We will need 
explicit formulae for the matrix elements /if^ of -f^2r^. The diagonalization of ifgr^ above 
gives one kind of formula. We also consider a second approach to inverting (due to 
[K]) via finite difference equations. The two approaches give quite different formulae for 
the inverse Hessian sums and have different applications in the inverse results. In several 
of the calculations in this section, we assume for simplicity of exposition that 7 is elliptic; 
the hyperbolic case is easier and all formulae analytically continue from the elliptic to the 
hyperbohc cases. 

For our purposes it will suffice to know the formulae for the elements h}^^.. To emphasize 
the fact that the matrix elements depend on, and only on, (r, a) we denote them by h^^.{a). 
The first formula comes directly from the diagonalization above. 
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where the denominators are defined in ([y3]j-(fy7|j. 

The second, finite difference, approach expresses the inverse Hessian matrix elements 
in terms of Chebychev polynomials T„, resp. ?7„, of the first, resp. second, kind. They are 
defined by: 

sin(n + 1)^ 



T„(cos^) = cos n9, f/„(cos^ 
Proposition 6.3. [KJ (p. 190) With the above notation, 



sin 9 



{-L)-'hlt{a) = 2[i,tJ(W2)] [U2r-,+p-ii-a/2) + [/,_p_i(-a/2)], 1 < p < q < 2r) 

We note that h^'^ = h'^^ so this formula determines all of the matrix elements. 

The special cases r = 1, 2 are already very helpful in the inverse problem. We recall that 

Ti{x) = X, T2{x) = 2x^ - 1, Tsix) = Ax^ - 3x, T^{x) = 8x^ - 8x^ + 1; 

Ui{x) = 2x, f/sO 
from which we calculate: 



(79) 
and 



(80) 



2x2- 


- 1, T3(x) = 


4x3 - 


3x, T4^{x] 


= 8x^ - 


= 4x2 


- 1, U;ix) -- 


= 8x3 


- 4x, U4 


= 16x2 - 




""^ 1 


f a 








a2-4 i 




a I 


V 






/a^ - 2a 


-a2 


2a 


-a2 \ 


L 


— a2 a^ 


- 2a 




-a2 


2a 


■4a2 


2a 


-a2 


a3 


- 2a 


-a2 




\ -a2 


2a 




-a2 


a3 - 2a) 



In terms of Floquet angles, we have (in the elliptic case). 



hence 



2[1 - Ta^ (cos a/2)] 



[t/2r-g+p-i(-cosa/2) + f/y_p_i (- COS a/2)] , (1 < p < g < 2r)), 



{-\)v-i ^sin(2-r-g+p)a/2 sin(g-p)a/2 -| 



(82) (-L)-i /i^" 



2[1— cos ra)] I- sinQ:/2 



sin Q:/2 



(1 < p < g < 2r) 



2r 



r sin(2r-p+g)a/2 , sin(p-g)a/2 1 
2[1— cosrajL sin«/2 sina/2 J 



'\<q<p<2r) 

We note that the expression in Proposition (16.21) is the Fourier inversion formula for ([82 
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COROLLARY 6.4. We have: (-L)-^ = ^l^ty, = = 2^ ^^t f . 

6.5. Linear sums. We now complete the proof of Corollary 15.111 by summing the matrix 
elements in the first row [H^^\ = (/ill, _ _ _ ^ f^ii2r)^ column) of the inverse. As a check on 
the notation and assumptions, we calculate it in two different ways: 

Proposition 6.5. Suppose that j is a 'L2-symmetric bouncing hall orbit. Then, for any p, 

^p2r ipq _ -L_ _ -L 
Z^g=l"'2r a+2 2-2coso/2' 

Proof. Because H2r^ is a circulant matrix, the column sum is the same for all columns. Hence 
we only need to consider the first column. 

(i) By Proposition 16. 2[ we have 

Y^2r , pq _ Y^2r _ ^ V^2r ■^2r-l to('?-l)fc 

Z^q=l "'2r — l^q=l 'hr — 2r Z^g=l Z^fc=0 pZAw^ 

= (—T] V^**"^ 4o ^ -L ^ ^ ^ -L 

\ ^) Z^k=0 pa,riw'') Pa.r{l) 2+a 2-2cosa/2' 

(ii) Since X^gLi ^It — Yl^Li ^2r ^^i' P = 1, • • • , 2r, we can set p = 1 in the sum over q to 
obtain, 

2r 2r 2r 

(83) 1= ^ v/^^^ = Ev]E^'']- 

p,q=l p=l g=l 

It then follows from (fT6|) and Proposition (12. 2p that (— L)^^ Ylp=i ^pq' — 2+0 = 2—2 cos a/2. 

□ 

6.6. Decoupling Balian-Bloch invariants. Corollary (15. lip expresses Byj-i in terms of 
inverse Hessian matrix elements. To prove Theorem II. 1^ it is essential to show that we can 
separately determine the two terms 

(1) (/^^K«))H2K^^.,o)/(2^)(0) +4^|^/(3)(0)/(2.-i)(0)}, 

(2) 4(n;.2. .3,0) j:lUhltia)f}f'^'KO)f'-'^-'KO)- 

2,j + l ^ 

To decouple the terms we prove that they have behave independently under iterates r of 
the bouncing ball orbit. We use the simple observation: 

Lemma 6.6. Let F3(r, a) = YllLii^lti^))^ ■ (^2r('^))~^-^3('^' ^■^ non-constant in r = 
1, 2, 3, ... , then both terms (l)-(2) can be determined from their sum as r ranges over N. 

Proof. Put 

(W(j2j-l,3,o) 

A = {2(^g2„o)/(2^)(0) + 4^2;^/(3)(0)/(2^-i)(0)}, B = 4(ti;^2.-M.o)/(3)(0)/(2^-i)(0). 

It is assumed that we know {h}^{a)YA + F3{r,a)B for all r G N. To determine A,B it is 
clearly sufficient that the matrix 

{hllia))' F3(r,a) 
{hllia))' F,{s,a) 
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is invertible for some integers r ^ s. But this says precisely that {h^{a))~'^F^{r^a) 7^ 
[h^g^a))'"^ F^{s ^ a) for some integers r ^ s. □ 

6.7. Cubic Hessian sums. We now prove that [h^^^a))''^ F^{r , a) is indeed non-constant 
for all but finitely many a. 

Proposition 6.7. The 'had' set B of [W\) consists of {0, -1, ±2}. 

Proof. We will give two different proofs of the finiteness of B. In both, we consider the sets 
Br,s = {a G M : {h^^,.{a))'^ F:,{r , a) = (/i^^ (a)) -2^3(3, a)}. 

6.7.1. First proof of Proposition \67^ Dedekind sums. The first is based on an explicit cal- 
culation of Fs{r,a) as a Dedekind sum. It is not very efficient in bounding the cardinality 
of Br^s but gives a clear proof that this set is finite. 

Lemma 6.8. We have : 

F(ra) = V - 

(2^)' ,^|;lo(a + 2cos^)(a + 2cos^)(a-t-2cos(^i±^)' 

In the hyperbolic case, we obtain a similar result with cos replaced by cosh . 

Proof. Using Proposition 16.21 we have (with w = , and = equal to congruence modulo 
2r), 

1 



0<fc,<2r-l;fci+fc2+fe3=0 (a+2cos !^){a+2cos ^)(a+2cos ^) 

= 2rV^''"^ i 

^ki,k2=0 ^^_^_2 cos ^1^ ) (a+2 COS ^ ) (a+2 cos ^M+MZ ) " 

□ 

We now complete the proof of Proposition 16.71 By Corollary 16.41 {hll{a))~^F3{r,a) is 

the rational function ( 2{i-t (-s^)) ) ^si^y^^)^ where as above, Tn,Un are the Chebychev 
polynomials. 

We now observe that for r ^ s, i^ 2{\-T {-^)) ) -^3(^5°) ^^"^ ( 2('i-~t\-^-'')) ) -^3('^''^) 
independent rational functions. Indeed, the poles for given r are the values a = — 2 cos | 
where ot = ^ for some = 1, . . . , 2r. Hence, there can exist only finitely many solutions of 
the equation 

(85) Ui-r.(-i)) j ^^''■'■'-i 2(i-r.(-|)) j ^^f^-"' 
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for any r ^ i.e. Br,s is finite. 

□ 

It is interesting to observe that the sums above are generahzed Dedekind sum, i.e. the 
sum ^^g£,^ -^3(C! A of ^^^6 function 



{z + cosxi)(2; + cosx2)(-2 + cos(xi + X2)) 



over the set of 2rth roots of unity ^mod 27rZ^ with k= {ki, /C2) G [0, 2r - 1] x [0, 2r - 1] 
of the torus. The summand is is a continuous periodic function of (a;i,X2) G [0, 1] x [0, 1] 
for z ^ [—1, 1]. In fact, I^{x^z) is also symmetric under inversion and refiection across the 
diagonal and the sum has additionally the form of a multiple Dedekind sum 

S2(l, 1; 2r) = /(^i' r)/(A;i + k^, r) 

k\ ,fc2(mod 2r) 

of two variables in the sense of L. Carlitz |Ca], with f{k,r) = (^^cos^Trfc/r) • 

We remark that under the non- degeneracy assumption that a/n ^ Q, cos a/2 is never a 
pole of F^lr, z) for any r. In the hyperbolic case, it is obvious that cosh a is never a pole of 
Fz{r,z). 



6.7.2. Second proof: Explicit inversion of the Hessian. We now give a second (and quite 
elementary) method of determining B by simply using the formulae for ( I79|) and H^^ 
(IHOj) . This calculation is due to the referee and to H. Hezari. 
From the explicit formula for H^^ we have: 



3 

3 



q=l 

Further, h}^ = -^zj- From the explicit formula for H^^ we have 

4 / r \ 3 



q=l ^ ' 

Further, hi' = J^^. 

Thus, Bi^2 is the set of solutions a of the equation 

a3_g (a2_4)2 _ (a4-4a2)2 a9-6a7-2a6+12a5 
(a2-4)3 a2 ~ (a3-2a)2 (a4-4a2)3 

(a^ - 2a)2(a3 - 8) = - 6a^ - 2a^ + 12a^ 

A little bit of cancellation reduces the equation to degree 6. The distinct roots are {0, —1, 2, —2}. 
QED 
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6.8. Final step in proof of Theorem ll.lt Inductive determination of Taylor coeffi- 
cients. We now prove by induction that on j that /^-' (O), f^^-^^^\0) are wave trace invariants, 
hence spectral invariants of the Laplacian among domains in Vi ^. 

It is clear for j = 1 since (1 — Lf^'^\0) = cos a/2 (resp. cosh a/2) and a is a Balian-Bloch 
(wave trace) invariant at 7 (see [^). In the case j = 2, the Balian-Bloch invariants have the 
form (Q. Using that a is a Balian-Bloch invariant and the decoupling argument of Lemma 
EDand Proposition [621 (/^^^O))^ is a spectral invariant. By reflecting the domain across the 
bouncing ball axis if necessary, we may assume with no loss of generality that f^^\0) > 0, and 
we have then determined {f^^\0)) from the sequence of Balian-Bloch invariants. Using again 
that a is determined by the Balian-Bloch invariants, it follows that /'•^-'(O) is determined. 

We now carry forward the argument by induction. As j J + l, we may assume 
that i7^-'~^/(0) is known. The terms denoted -R2ri7^"'~^/(0) in Theorem 15.11 are univer- 
sal polynomials in the data i7^"'~^/(0), hence are also known. Thus, it suflices to determine 
/(^■'^(O), f^'^^~^\0) from (jni). By the decoupling argument, we can determine {f^^\0)){f^'^^~^^{0)), 
hence {f^^^-^\0)), as long as {f^^\0)) ^ 0. But then we can determine f^^-'\0). By induction, 
/ is determined and hence the domain. 

This completes the proof of Theorem ( 11. ip . QED 

Remark 6.9. From this argument it is only necessary that the coefficients wg etc. are 
non-zero and universal. It is not necessary to know the precise values of the coefficients of 
/(2^)(0),/(2^-i)(0). 

6.9. The case where f^^^O) = 0. If f^^\0) = 0, the inductive argument clearly breaks 
down. There is a natural analogue of it as long as /'•^•'(O) 7^ 0. We only sketch the analogue 
to make it seem plausible, but do not provide a complete proof. 

Instead of inductively determining f^'^^\0), f^'^-'~^\0), we inductively determine f^^^\0), f^'^-'~^\0) 
by a similar argument. Since /(3)(o) = 0, the terms f^^^'^^O) have zero coeflicients, and 
each new 'odd' term as j — j + 1 now has the form [X]g=i(^^^)^]/''^HO)/''^"'~'^''(0)- To carry 
out the analogue of the previous argument, it suflices to show that ^2r^[X]g=i(^^^)^] is a 
non-constant function of r. It should be plausible that this is the case, at least if we exclude 
a flnite number of values of the Floquet exponents. 

There then arises an inflnite sequence of further sub-cases where all odd derivatives vanish 
up to some jo + 1. To handle this case, we would need to show that /i2r^[X]g=i(^^^)^"'"'''"^] is 
non-constant for all jo- This should again be plausible. 

In the case where all odd derivatives vanish, the function /_|_ is even and the proof reduces 
to the previously established case of two symmetries. 

7. Proof of Theorem (11.41) 

We now generalize the results from a bouncing ball orbit to iterates of a primitive Dm- 
invariant m-link reflecting ray 7. For short, we call 7 a Dm-ray. 

7.1. Structure of coeflflcients at a Z^m-ray. 

7.1.1. Dm-rays. In the dihedral case, we orient Vl so that the center of the dihedral action 
is (0, 0) and so that one vertex vq of 7 lies on the y-axis. We again deflne a small strip 
Tf;(7), which intersects the boundary in n arcs. We label the one through vq by a. We then 
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write a as the graph y = f{x) of a function defined on a small interval around (0, 0) on 
the horizontal axis. Since we are only considering Z)„-invariant rays, the domain is entirely 
determined by a and /. 

We first need to choose a convenient parametrization of dQ fl T^{'y)- Either a polar 
parametrization or a Cartesian parametrization would do. For ease of comparison to the 
bouncing ball case, we prefer the Cartesian one. Thus, we use the parametrization x e 
(— e, e) — > (x, f{x)) for the a piece. We then use x — > -R27r/m(-^' /(^)) rotate Rl-^/jj^a. 

When considering Y, we need variables Xjs{j = 1, . . . , m; s = 1, . . . , r),Xjs — > Riw/mi-'^js, xjs)- 
We have: 

R;^:/jx,j{x,)) = {xp'\{f{x,)r^^)) 

:= (cos(2p7r/m)xp + sm{2pn / m) f (xp) , — sin(2p7r/m)xp + cos(2p7r/m)/(xp)). 

We also put {-lj'{x,)r(p) .= Rll^/j-l, f'{x,)). 
We then define the length functional 

C''{y,Xo,Xi,. . . ,Xmr) = \{Xo,y) - (Xi, f {Xi)Y^'^^ \ + |(Xo,l/) - (Xrm, f (Xmi))''^''"'^ \ 

(86) 

\{x,j{x,)r^^^ - {x,^uf{x,^i)r^^^'^\ 

We will need a formula for its Hessian in the case of a Dm-ray- By ( |KTj . Proposition 3), 
the Hessian Hrm in x — y coordinates at the critical point (xi, . . . , Xrm) corresponding to Y 
is given by the matrix (\2.2\j with s = -^^j- 

A key point in what follows (as in |Zlt IZ2j ) is that the reflection symmetry of a and / 
implies that f^'^^~^\0) = for all j. This eliminates the most serious obstacle to recovering 
/ from the wave trace invariants at 7'', namely the fact that in the transition from the jth 
Balian-Bloch invariant to the (j + l)st, two new derivatives of / appear. 

As in the Z2-symmetric case, there are principal and non-principal terms. The principal 
term in the Dm case, analogously to the bouncing ball case, equals Trp * N^"' o N[{k) o x{k) 
for r repetitions of the dihedrally symmetric orbit. 

In analogy to Lemma (15.11) we prove: 

Lemma 7.1. Let 'y be a Dm- ray, and let p be a smooth cutoff to t = rL^ as above. Then: 

• B^rj = Pm,r,jU'^'^K^), /^^^O), ■ ■ ■ , Z^^^'+^HO)) '^here P2,r,i(^i, • • • ,6i) «s a polynomial. 
It is homogeneous of degree —j under the dilation f A/, is invariant under the 
substitution f{x) f{—x), and has degree j + I in the Floquet data e*"*". 

• In the expansion in Theorem (1.1) of |Z5j of TrRp{{k + ir)) , f^^^\0) appears first 
in the k~^^^st order term, and then only in the k~^^^st order term in the stationary 
phase expansion of the principal term Trp * N^^^ o N[{k) o x{k); 

• This coefficient has the form 

B,rj_, = mr{h'y f(^^\0) + Rmr{J''-'fm, 
where the remainder RmriJ^^-'"^ f{0)) is a polynomial in the designated jet of f. 
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We use the analogue of Theorem 13.11 for the case of the dihedral ray. As in the case 
of a bouncing ball orbit, we have a finite number of oscillatory integrals arising from 
the regularization of the trace. We express the resulting oscillatory integrals in Cartesian 
coordinates of (polar coordinates are also convenient for this calculation). We put x = 
{xo,yo). Each oscillatory integral J^}^ localizes at critical points, we may insert a cutoff to 
Te(7). This gives m*^ possible terms, corresponding to the possible choices of the arcs in the 
product {dn n T,(7))*^. We put: 

{m^'}:={a:ZM^{l,...,m}}, 

and write 

:= (cos(2p7r/m)xp + sin(2p7r/m)/(xp), — sin(2p7r/m)xp + cos(2p7r/m)/(a;p)). 

We also put {-lj'{x,)r(p) .= r;^J/J-1J'{x,)). 

The oscillatory integrals have the phase functions on [dQ fl T^{'j)Y'^ of the form: 

mr—l 

(87) C''{x,,...,Xmr) = \{^pJi^p)r^'^ - {Xp+,,f{Xp+,) f ^+'^\ 

Only 2m cr's ( 2 modulo cyclic permutations) give length functions which have critical 
points with critical value rL^, namely the ones ctq where (Jo{n) = R{±n27T /m). Indeed, 
the only Snell polygon with this length is Y by assumption, and so {x1^^\ . . . ,Xrm"^^) must 
correspond to the vertices of 7^*". Since the good length functions represent isometric situa- 
tions, it suffices to consider the case ao{n) = R{n2Tr/m). In this case, we denote the length 
function simply by L and to simplify the notation we drop the subscript in cxo. 

We now make a stationary phase analysis as in the bouncing ball case to obtain the 
expressions in Theorem (15. ip . As mentioned above, there are two principal terms: The 
principal oscillatory integrals I^m^p^ those in which M = rm and in which no factors of 
Nq occur, i.e. cto(j) = 1 for all j = 1, . . . ,rm. Also, there are now m components of the 
boundary at the reflection points, and w± cycles around them for r iterates. 

7.2. The principal terms. They have the phase 

mr—1 



and the amphtude 

(89) 



2 



a^k, Xi, . . . , Xrnr, y) = U;LMik + it) ^(x^i^'^ " (^^P^"^)' + {f {x,-,Y^P-^) - /(x,)'^^)^) 
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We observe that it has the form A{x, y, f, /'). The /' dependence will be particularly impor- 
tant later on. 

7.2.1. The principal term: The data /^•'(O). As in the bouncing ball case, by the same 
argument, the data /'■^-'^(O) appears first in the term of order k~^^^ and it appears linearly 
in the term a^WR^. We now show that its coefficient is given by the formula in Lemma 
(17. ip . Due to symmetry, it suffices to consider any axis and one endpoint of it. We observe 
that only the 'diagonal terms' of W , i.e. those involving only derivatives in a single variable 
can produce the factor /'•^■'-'(O). Since /'(O) = x\x=q = and since the angle between 
successive links and the normal equals vr/m an examination of (H5ll shows that the coefficient 
of /(2i)(0) equals 

The data f^'^-'^^\0) vanishes due to the symmetry around each dihedral axis. 

Finally, as in the bouncing ball case, and for the same reasons, non-principal oscillatory 
integrals do not contribute to this data. 

This completes the proof of Lemma (17. ip . □ 

Remark It would also be natural to employ polar coordinates in the proof. In that case, we 
align fl so that one of the reflection axes is the positive x-axis, and express dfl parametrically 
in the form r = r{'&) where is the angle to the x-axis. Then r(—'&) = r(-t}),r(-t}+^) = r (■(}). 
The goal then is to determine r. To do so, we write out that q{'d) = (r('i?) cos('i9), r(^9) sin('(9) 
and compute as above. We flnd that r^'^^\0) arises flrst in the k~^^^ term with the same 
coefficient as for /'^^■'•'(O) above. The rest of the proof proceeds as with Cartesian coordinates. 

7.3. Dihedral domains: Proof of Theorem ( II. 4h . We now complete the proof of The- 
orem (11.41) . 

We prove by induction on j that /^■'(O) is a Balian-Bloch invariant. It is clear for j = 1 
since (1 — Lf^'^\0) = cos{h)a/2 and a is a Balian-Bloch (wave trace) invariant at 7. In 
general, the eigenvalues of are wave trace invariants [F]. 

Assuming the result for ri < j — 1, it follows that p^n*-! is a spectral invariant. It thus 
suffices to extract p^{0) frompO^-.^, i.e. from {EpLi(^^^)^}/^^^HO)- 

Thus, the only missing step is to show that if 7 is Dm-ray, then the h'^^ are Balian-Bloch 
invariants of 7*". In other words, that s is a wave trace invariant. If A, denote the 
eigenvalues of P^r, then we have A + A^^ = 2 + det Hmr- Here we use that all bj equal 1. It 
follows that s is a function of A, hence that it is a Balian-Bloch invariant. 

The proof of Theorem (11.40 is complete. 
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